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Preparation of D-Arabinose-l-C* and D-Ribose-1-C’ 


Harriet L. Frush and Horace S. Isbell 
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b-ribonate-1-C", 


p-arabinose-1-C™ and p 
of 30 and 8.5 percent, 


b-erythrose 


radiochemical vields 


General acid catalysts in the cyanohydrin reaction appear to favor formation 


resulting from the reaction of labeled cvanice 
were separated 
respectively 


crystalline 
rhe were 


as potassium bD- 


salts converted 


to the corresponding lactones, and these were reduced to the sugars by use of sodium amalgam 


in the presence of sodium acid oxalate 


1. Introduction 


\s part of a program to make position-labeled 
; available for research workers in other labora- 
methods have been developed at the National 
. of Standards for the preparation of p-glucose- 
p-mannose-1-C', p-mannitol-1-C™, p-fruc- 
6-C™, lactose-1-C™, and p-arabinose-5-C"™ [1, 
» 3, 4, 5, 6)2 The present report gives methods 
for the preparation of p-arabinose-1-C™ and p-ribose- 
1-{ p-Arabinose-1-C' was required as an inter- 
mediate in the preparation of p-glucose-2-C™, and 
both p-arabinose-1-C" and p-ribose-1-C" were needed 
for st of the transformation of pentoses in 
biological systems. Rappoport and Hassid [7] 
obtained t-arabinose-1-C'™ by application ‘of the 
Sowden-Fischer nitromethane synthesis [8] to L-ery- 
throse. The t-arabinose-1-C was separated in 
3-percent vield by partition chromatography 
In light of prior work at the Bureau, it seemed de- 
sirable to attempt the synthesis of both p-arabinose- 
1-C" and p-ribose-1-C™ by the cyanohydrin method, 
beginning with p-erythrose. It was found that in 
the reaction of cyanide with p-erythrose the presence 
of a general acid catalyst, such as bicarbonate or 
ammonium ion, favors formation of the arabonic 
epmer. The epimeric products of reaction were 
separated by the following steps: (1) Hydrolysis of 
the nitriles with aqueous sodium carbonate, (2 
passage of the solution over a cation exchange resin 
the free acids, (3) neutralization of the acids 
with potassium hydroxide, and separation of labeled 
mic acid as crystalline potassium p-arabonate, 
version of the potassium salts in the mother 
to cadmium salts by passage of the solution 
cation exchange resin, and neutralization of 
the acid with cadmium hydroxide, and (5) separation 
of the labeled p-ribonic acid in the form of crystalline 
cadmium p-ribonate. 
This procedure gave labeled potassium p-arabonate 
and cadmium p-ribonate in radiochemical yields of 
and 23 percent, respectively. The potassium 
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p-arabonate-1-C'™ was converted to p-arabono-y- 
lactone-1-C™, and this was reduced with sodium 
amalgam in the presence of sodium acid oxalate 
A 56-percent vield of the crystalline p-arabinose-1-C" 
was separated without carrier; by use of carrier, the 
radiochemical yield was increased to 60 percent 
As potassium p-arabonate-1-C™ was produced in 50- 
percent yield, the over-all radiochemical vield of 
p-arabinose-1-C' was 30 percent, based on the 
evanide originally used 

The cadmium p-ribonate-1-C"', prepared from the 
mother liquors of the potassium p-arabonate-1-C", 
was converted to p-ribono-y-lactone-1-C™, which was 
brought to crystallization by nucleation. The 
lactone was reduced with sodium amalgam and gave 
p-ribose-1-C™ in 37-percent vield. Inasmuch as 23 
percent of the activity of the cyanide had been 
obtained as cadmium p-ribonate-1-C™, the over-all 
radiochemical yield of p-ribose-1-C' was 8.5 percent 
Work still in progress will undoubtedly raise the 
vields of both cadmium p-ribonate and p-ribose-1-C"™ 


2. Experimental Procedures 


2.1 Preparation of D-Erythrose 


For the preparation of sirupy b-erythrose, Sow- 
den’s method of periodate oxidation of a 4,6-substi- 
tuted glucose [9] was applied to 4,6-ethylidene-p- 
glucose essentially as described by Rappoport and 
Hassid [7], but with certain convenient modifications 
A mixture of 5.64 g of 4,6-ethvlidene-p-glucose, mp 
180° to 181° C [10] and 4.6 g of sodium bicarbonate 
was dissolved in an ice-cold solution containing 11.7 g 
of sodium metaperiodate in 150 ml of water. The 
solution was kept in an ice bath for a few minutes, 
and finally at room temperature for 2.5 hr; it was then 
freeze-dried The fluffy residue was extracted with 
a total of 150 ml of hot ethyl acetate in 3 portions. 
The extract was filtered through a bed of decolorizing 
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carbon and diatomaceous earth, and the solvent was 
removed by distillation under reduced pressure. 
The residue, amorphous 2,4-ethylidene-p-erythrose, 
was hvdrolyzed by refluxing for 1 hr with 60 ml of 
0.1-N sulfuric acid; the hydroly zate was cooled, and 
deionized by passage through a column of mixed 
eation ® and anion ® exchange resins. The resulting 
solution was concentrated under reduced pressure to 
about 50 ml, and was then freeze-dried. The color- 
less sirup weighed 2.91 g and had a specific rotation, 
al”, of —17.3 The cyanide-combining power of 
the sirup was about 90 percent of the theoretical, 
and consequently it was assumed that the residue 


contained 90 percent of p-erythrose It was used 


without further purification for the cyanohydrin 
svnthesis described in the following sections 


2.2. Effect of General Acid Catalysts in the Cyan- 
ohydrin Synthesis on the Yield of the D-Arabonic 


Epimer 


Several cyanohydrin syntheses employing C"- 
labeled cyanide and sirupy pv-erythrose were carried 
out on a semimicro scale under the conditions de- 
scribed below. After hydrolysis of the cyanohydrins, 
each mixture was treated with a large excess of non- 
radioactive potassium p-arabonate. From the ra- 
dioactivity of the recrystallized potassium salt, the 
proportion of p-arabonic acid was calculated. 

The mixtures, two of which are listed in table 1 as 
experiments | and 2, were frozen, and the tubes’were 
sealed and allowed to stand at room temperature for 
72 hr. They were then opened and heated on a 
water bath at 80° C for 7 hr in the presence of a 
stream of air. A few drops of water were added 
from time to time, but at the end of the hydrolysis 
period the solutions were allowed to evaporate to 
dryness. The residue in each tube was dissolved in 
5 ml of water, and 500 mg of potassium p-arabonate 
was added. The mixture was warmed slightly to 
dissolve the salt, and the solution was filtered with 
the aid of a small amount of decolorizing carbon, and 
treated with 3 volumes of ethanol. Potassium p- 
arabonate crystallized freely from all preparations. 
After storage for 18 hr in a refrigerator, the mother 
liquor was removed from the crystals, which were 
washed in place with cold, aqueous ethanol (1:4). 
The samples were recrystallized once from water and 
ethanol, dried in a vacuum desiccator, and assayed. 
Experiments 1 and 2 reported in table 1 indicate 
that ammonium carbonate favors 
p-arabonic nitrile. A large-scale preparation re- 
ported in section 2.3, and also listed as experiment 
3 in table 1, shows that sodium bicarbonate likewise 
favors production of the arabonic epimer. The 
results are in accord with the prior observation that 
general acid catalysts alter the proportion of the 
epimeric nitriles [3, 5). 
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TABLE |. 
arabonic 
D-¢ rythrose 


Effect of a general acid catalyst on the yield 


epimer obtained by the addition of cyar 


Activity 
found in 
carrier 
potassium 
p-arabon- 
ite > 


Reaction mixture * 


Millimoles 

°).0267 p-erythrose 

0134 NagCOs 
¢.0267 NaCN 
¢.0267 D-erythrose 

134 (NH4)sCOs 
¢.0267 NaCN 
4.0 p-erythrose 


* Total volume was 1 ml in experiments 1 and 2, and 2) ml in expe 
» Activities determined by means of a vibrating-reed electromet« 
modified Van Slyke-Folch wet oxidation procedure (see [11] 

Assuming 90-percent purity of the erythrose sirup 

* Calculated for entire 500 mg of carrier added 

¢ Contained 26.7 wc of carbon 14 

t Radioactivity recovered in several crops 

« Contained 14 me of carbon 14 


2.3. Preparation of Potassium D-Arabonate-l-C" fron 
D-Erythrose by the Cyanohydrin Synthesis in th 
Presence of Sodium Bicarbonate 


Ten milliliters of an aqueous solution contain 
3.1 millimoles of sodium cyanide (14.0 me of carl 
14) and 5 millimoles of sodium hydroxide was froz 
on the sides of a small glass-stoppered tube 
small lump of solid carbon dioxide was added, a 
then a solution containing 5 millimoles of sodiu 
bicarbonate and approximately 4 millimoles of 
erythrose in i0 ml of water. The loosely stopper 
flask was placed in an ice bath, shaken until t! 
contents had dissolved, and then kept in the ba 
for 1 day, and at room temperature for 5 days 
The nitriles were hydrolyzed by heating on a stea 
bath for 7 hr, and the resulting ammonia ws 
removed by evaporation of the solution to drynes 
in a stream of air. The residue was taken up 
water and passed over a column containing 10 ml ¢ 
cation exchange resin’? to remove sodium ions. Tl 
column was washed until the activity of the effluer 
was negligible, and the effluent was concentrat 
under reduced pressure to about 25 ml and neutral. 
ized with potassium hydroxide. Two grams 
nonradioactive potassium p-arabonate was dissolv: 
in the sirup and evaporation was continued to ! 
volume of about 4 ml; the solution was then satv- 
rated with methanol. After 18 hr, a crop of potassiu: 
p-arabonate-1-C" was obtained, which, when rr 
crystallized and dried, weighed 2.23 g and had 4 
activity of 5,980 we. By the further addition 
2.5 g of carrier potassium p-arabonate in 3 portions 
970 ye of activity was obtained. Hence the told 
radiochemical yield of potassium p-arabonate-1l- 
was 6,950 uc, or 50 percent of the activity of th 
original cyanide. The residual sirup was used {« 
the preparation of cadmium p-ribonate-1-C" (se 
sec. 2.6). 

7 Amberlite IR-120, analytical grade, Resinous Products Division of hoht 
Haas Co., Philadelphia, Pa 
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Preparation of Potassium D-Arabonate-l-C" from 
Erythrose by the Cyanohydrin Synthesis in the 
»sence of Ammonium Carbonate 


another p-arabinose-1-C™, the 
nium carbonate method was emploved. Five 
ters of a solution containing 0.7 millimole of 
m hydroxide and about 0.7 millimole of sodium 
le with 5.4 me of carbon 14 was frozen on the 
ff a small glass-stoppered tube. Five milli- 
of approximately 1 \f ammonium carbonate 
hen frozen in the tube, and finally 1 ml of a 
ym containing 1 millimole of p-erythrose*® was 

The mixture was allowed to thaw in an ice 
and was then kept at temperature for 
s. It treated the manner 
bed in section 2.3, except that a crop of potas- 
p-arabonate-1-C™ (52 meg, containing 2,030 
vas separated without the addition of carrier. 
otal radiochemical vield (2,446 ue, or 45 per- 
was slightly than that of the former 
ation. However, as greater mechanical 
curred in this preparation because of the higher 
of activity, it is believed that the two methods 
qually satisfactory in the synthesis of potas- 
p-arabonate-1-C }*, 


preparation ol 


room 


was essentially in 


less 


loss 


Preparation of D-Arabono-y-Lactone-l-C™ and 
Reduction to D-Arabinose-l-C" 


solution containing 446 mg (2.19 
of potassium p-arabonate with-1,190 ye 
tivity was passed through a column of cation 
exchange The column was 
iwhly washed, and the effluent was evaporated 
0° C substantially to dryness. The residue, 
olved in a few milliliters of methanol, was trans- 
d to a reduction tube such as that described 
ously [3]. Briefly, it consists of a test tube 2.8 
20 em fitted with a stirrer reaching to the bottom 
tube and with a sidearm for the introduction 
amalgam. The solution in the tube was con- 
ated in an air stream and seeded with p-arabono- 
tone."° The tube was finally stored in a desic- 
and the contents was moistened from time to 
with methanol. When lactonization was com- 
as indicated by the disappearance of any 
sirupy phase (1 or 2 weeks), the material was re- 
| with sodium amalgam by the procedure 
cribed later. 
lhe laectonization described above required con- 
siderable time. In order to avoid this delay, the 
following rapid procedure was developed: Two 
millimoles of potassium p-arabonate-1-C' was passed 


aqueous 


noles 


resin (see footnote 7 


of uncertainty in the concentration of the cyanide solution, the 
was in larger excess than would ordinarily be added 
t for the analyses in table 1, all assays of carbon 14 were made in form- 
ition [12]. Details of the analysis of this crop of potassium D-arabonate- 
llustrative: The material was dissolved in water and the volume was 
10 ml A 50-\ aliquot was diluted to 1 ml with water in a micromixing 
A 20-\ aliquot of this was diluted to 1 m! with formamide, and the 
gave 74.9 counts per second in a 2-* proportional counter in which 1 
second is equivalent to 0.00271 we 
ree acid may crystallize from the freshly prepared sirup at this point. 
1 of the acid to the lactone at room temperature requires at least a week, 


through a column of cati yn-exchange resin, and the 
effluent was concentrated under reduced pressure to 
a sirup, which was transferred with a few drops of 
water to a reduction tube. Five milliliters of glacial 
acetic acid was added, and the tube was placed in a 
boiling-water bath and heated in a stream of air 
for 5 hr. After 1 hr of this period, 1 ml of water and 
5 ml of were added \t three 1-hr in- 
tervals thereafter, the material was treated with 5 ml 
of acetic acid only \t the end of the heating period, 
the pale-vellow residue was dissolved in 3 or 4 drops 
of methanol, seeded with p-arabono-y-lactone, and 
placed in a desiccator over calcium chloride. It 
crystallized well within 24 hr 

When samples of nonradioactive p-arabonic acid 
were subjected to this procedure, the optical rota- 
tions, [a]?’, 24 and 72 hr after the sirups were seeded 
with the crystalline were, respectively, 
-66.0° and +67.6 the weight of the 
lactone These values correspond approximately 
to conversion of 91 and 94 percent, respectively, of 
the p-arabonic acid to the lactone 

Prior to the preparation of the labeled sugar, a 
sodium amalgam reductions was carried 
ut on the nonradioactive lactone in order to de- 
termine optimum conditions. Sodium acid oxalate 
was employed as a neutralizing agent [6], and after 
removal of the salts in the manner described below, 
the reduction product was analyzed for reducing 
sugars by the modified Scales method [13]. It was 
found that a maximum yield of sugar (85 percent) 


acetic acid 


lactone 
based on 


series of 


is obtained with 4.6 g of 5-percent sodium amalgam 


per millimole of lactone 

Ina ty pical reduction of the labeled lactone to the 
sugar, the reduction tube, containing 6.4 g of crystal- 
line sodium acid oxalate and the lactone prepared 
from 446 mg of potassium D-arabonate-1-C ", having 
1,190 we of activity, was placed in an ice bath. After 
the addition of 20 ml of ice water, the stirrer was adjust- 
ed to produce vigorous stirring at the bottom of the 
tube, and 9.2 g of 5-percent sodium amalgam pellets ™ 
was added through the side arm. Stirring was con- 
tinued for 3 hr, at the end of which time the mercury 
was removed, and the aqueous mixture was diluted 
with 5 volumes of methanol. The crystalline salts 
were separated, washed with methanol, and dis- 
carded, after it was ascertained that they were in- 
active. Phenolphthalein indicator was added to the 
mother liquor, previously cooled in ice water, and 
then aqueous sodium hydroxide until a faint but 
permanent pink color was obtained. The solution 
was concentrated under reduced pressure to about 10 
ml, and treated with 5 volumes of methanol; the pre- 
cipitated salts were separated, washed with methanol, 
tested for radioactivity, and discarded. The solu- 
tion was concentrated to remove the alcohol, and 
passed over a column: containing 25 ml of mixed 
anion and cation exchange resins (see footnotes 5 and 
6); the resins were washed until the activity of the 
washings was negligible. The effluent was con- 

The amalgam was prepared in small pellets by pouring it in molten condition 
through a heated alundum thimble, having a small hole in the bottom, into a 


2-ft “shot tower” of mineral oil The pellets were stored under mineral oil, and 
just before use, were blotted dry, weighed, and rinsed with benzene 
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centrated under reduced pressure and at a tempera- 
ture less than 40° C to about 10 ml; at this stage the 
solution was filtered into a round-bottomed flask 
through a microfilter containing a small bed of de- 
colorizing carbon,” and the combined solution and 
washings were freeze-dried in a 50-ml flask. The 
residue in the flask was dissolved in a few drops of 
methanol, and 2-propanol was added to incipient 
Crystallization began immediately after 
At the end of 2 days, the mother liquor 


turbidity 

nucleation 
was separated, and the crystals were washed in place 
with a few drops of a methanol-2-propanol mixture 
2:1). For reerystallization, the crude material, dis- 


solved in a few drops of water, was transferred to a 
standard-taper test tube, and the water was removed 


under reduced pressure. The thick sirup was dis- 
solved in a few drops of methanol, and 2-propanol 
was cautiously added, until saturation was ap- 
proached. The crystals that formed in 24 hr 
weighed 185 mg and contained 672 ye of activity. 
With the aid of 100 mg of nonradioactive p-arabinose 
as carrier, an additional 39 ye of activity was ob- 
tained from the mother liquor. The yield of p-arab- 
inose-1-C™ (711 we) from potassium pD-arabonate-1- 
C" was thus 60 percent 


2.6. Separation of Cadmium D-Ribonate-1-C" 


The mother liquor from the preparation of potas- 
sium p-arabonate-1-C"™, described in section 2.3, was 
evaporated to remove the alcohol present. The 
residue was dissolved in water and passed over a 
column of cation exchange resin (see footnote 7). 
The effluent, containing 6,900 ye of activity, was 
then neutralized with cadmium hydroxide to the end 
point of phenolphthalein. After the addition of 1 ¢g 
of nonradioactive cadmium p-ribonate, the solution 
was filtered with the aid of decolorizing carbon, and 
the filtrate was evaporated under reduced pressure 
toa thin sirup. Nucleation of the sirup with erystal- 
line cadmium p-ribonate and careful addition of 
aqueous methanol resulted in a voluminous crystal- 
lization. The crystals were separated by filtration, 
and recrystallized by dissolving in 4 ml of hot“water 
and adding successively 5 ml of methanol and water 
(1:1) and 2 ml of methanol and water (2:1 The hot 
solution was then saturated with methanol and al- 
lowed to cool slowly. The first crop of recrystallized 
cadmium p-ribonate-1-C™ weighed 678.6 mg and had 
an activity of 1,796 we. After successive additions 
of two 1-g quantities of nonradioactive cadmium 
p-ribonate as carrier, and recrystallization of the 
resulting crops of salt, an additional quantity of 
cadmium p-ribonate-1-C™ containing 1,416 uc Was 
recovered. The total vield of cadmium p-ribonate-1- 
C™ was thus 3,212 ue, corresponding to 23 percent 
of the C™-labeled cyanide used 


This step rem 
amalgam 

‘In other 
obtained, 


preparations 


2.7. Preparation of D-Ribono-y-lactone-l1-C" and 
Reduction to D-Ribose-1-C'* 


A solution containing 679 mg of ecadn 
p-ribonate-1-C™, with an activity of 1,800 ye, 
passed over a column of cation exchange resin 
footnote 7), and the effluent was concentrated w 
reduced pressure to a sirup. The sirup was tr 
ferred by use of water to three reduction tubes 
evaporated to dryness by heating on a steam | 
under a stream of dry air. The residue in each | 
was dissolved in a few drops of water and the s 
tion, after the addition of 5 ml of methy! celloso 
was evaporated again. This process was repea 
4 times over a period of 2 days. Each residue w; 
then dissolved in 2 ml of methyl! cellosolve, and 
solution was seeded with crystalline p-ribono-y-lar- 
tone and stored in a desiccator over calcium chlorid; 
From time to time the residues were moistened wit! 
methyl cellosolve. Crystallization occurred slowl 
but ultimately all sirup disappeared, and the residues 
appeared to be entirely crystalline. 

A series of reductions on nonradioactive p-ribono 
y-lactone showed that a maximum yield of suga 
(70 percent) was obtained when 9.2 g of 5-percen 
amalgam and excess sodium acid oxalate were em- 
ployed for the reduction of 1 millimole of lacton 
Each of the three samples of lactone mentioned aboy 
was dissolved in 20 ml of ice water and reduced } 
use of the above quantities of amalgam and oxalat: 
The product from the three reductions was combine 
the mercury was separated, and the mixture was 
diluted with 5 volumes of methanol. The resulting 
crystalline salts were separated by filtration, washe 
with methanol, and discarded after a test show 
negligible radioactivity. The filtrate was conce! 
trated under reduced pressure to a volume of abo 
15 ml and diluted with 75 ml of methanol; the sodiun 
salts that separated were removed by filtration an: 
washed with methanol. As a test showed that this 
second crop of salts likewise contained no appreciab) 
radioactivity, it was discarded. The alcoholic filtrat 
was concentrated to about 10 ml, neutralized wit! 
aqueous sodium hydroxide in the presence of pheno 
phthalein until a permanent pink color was obtaine: 
and immediately deionized by passage over a colum: 
of mixed cation and anion exchange resins (see foot- 
notes 5 and 6). The effluent and wash liquor wer 
combined and concentrated under reduced pressu: 
to a volume of 6 ml; 750 mg of carrier p-ribose was 
dissolved in the sirup, and the solution was filter: 
through 1 ml of a mixture of decolorizing carbon an 
diatomaceous earth. The filtered solution was freez 
dried, and the residue was dissolved in 0.5 ml o 
ethanol. A mixture of ethanol and isoamyl alcoho 
was added to the point of incipient turbidity, and th 
solution was seeded with p-ribose and stored in : 
refrigerator. When crystal growth seemed complet 
(5 days), the mother liquor was removed wit! 
capillary pipette, and the crystals were carefull) 
washed, in place, with ethanol and dried over calewun 





768 mg) had an activity 


de The crude product 
7 wc, or 37 percent of the activity of the cadmium 
mate—1—-C"™ 

recrystallization, the product was dissolved in 
drops of water, a drop of acetic acid was added 


used 


through decolorizing 
under re- 
residue was 


filtered 
was 


the solution was 
The filtrate concentrated 
pressure to a thick sirup. The 
ved in ethanol, and the solution 
again under reduced pressure. The residue 
ssolved in 0.5 ml of ethanol, and after addition 
imvl alcohol to incipient turbidity, the solution 
eeded and stored, first at room temperature and 
in a refrigerator After 5 davs, the product 
eparated, washed in place with a few drops of 
ol and dried The recrystallized p-ribose-1-C 
i) contained 465 ue of carbon 14 By use of 
s on the mother liquor from the recrystalliza- 
radioactivity of the crude 


was concen- 
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Corrosion of Nickel Cast Irons in Soils 


Irving A. Denison and Melvin Romanoff 


Ihe 7 iits Oo ‘ sure! nts of the corros 
I ts conditions for 
orrosion, as det 
composition of t 
specimens were 


removed from tl OTTos\t\ 


» hydraulic press 


1. Introduction 
osion of cast iron in soils is characterized by 
velopment of an adherent laver of corrosion 
ts that increases In 
sses. The thickening of this laver is due in 

electrochemical action between the ferritic 
iphitic constituents of cast iron and partly to 
Sin pote ntial that arise from contact of th 
tal with the soil, as, for example, by differen- 
ration To the extent that corrosion of cast 
caused by differences in potential within the 
tself, improved corrosion resistance would bs 
d to result from improvement in the quality 
casting. Wesley, Copson, and LaQue [1] 
that the addition of small amounts of 
nd nickel plus chromium alters the structurs 
iron in such manner that calvanit action 
the different constituents is reduced 
lv. Larger additions of these alloying ele- 
n amounts sufficient to produce an austenitic 
shown to be considerably more effec- 
reducing corrosion under the experimental 
of the ennobling effect of nickel 
omium on the pote ntial of iron 
der to evaluate the effect of additions of nickel 
nickel plus chromium on the corrosion of cast 
soils, samples of alloy cast irons were included 


thickness as corrosion 


nown 


con- 


were 


ons he cause 


xtensive series of exposure tests of materials 
lerground construction. In 1932, samples of a 
lov austenitic cast-iron pipe were buried at 
t sites, and in 1941 samples of several low-alloy 
ons were buried at 13 of th 
| other site. In a previous report [2], corro- 
ita were reported for five periods of exposure 

samples of high-alloy cast iron and for a 
period for the samples of low-alloy cast irons. 
the completion of the exposure tests in 1952, 
rv the latter materials became available for 
riods of exposure, with a maximum exposure 


sites, 


same test 


ears 

removal of the samples from the test sites, 
wedure followed during the early p: riods of the 
ire program was to clean the samples of all cor- 
products preparatory to measuring the weight 
and the depths of the deepest pits. This pro- 


brackets refer to literature refer« 


» 5OO Ib 


cedure was later modified by subjecting those samples 
in the form of pipe to hydraulic pressure before re- 
moval of the corrosion products. This modification 
was introduced in order to evaluate the strength of 
the corrosion products remaining in the pits of the 
samples that had been perforated by corrosion. The 
methods used for cleaning the corroded samples and 
for measuring the depth of pitting have been pre- 
viously described [3] 


2. Properties of the Soils at the Test Sites 


The location of the t sites, the identification of 
the soil types, and the properties of the soils that are 
associated with corrosion are given in table 1. The 
retentiveness of the water is indicated by 
the values for moisture equivalent, which is the 
quantity of water retained by a previously saturated 
soil against a centrifugal force of 1,000 times gravity 
Because the real specific gravities of rock-forming 
minerals lie within a narrow range, values for ap- 
parent specific gravity are an index of the relative 
compactness or porosity of inorganic soils 


soils for 


3. Description of the Materials 


The forms, dimensions, and compositions of the 
specimens are given in table 2. The ends of the 
pipes were closed with cast-iron screw caps of the 
same composition as the specimens in order to exclude 
moisture from the interior. As an extra precaution 
against internal corrosion, the interior surface of the 
specimens was coated with heavy grease. 


4. Weight Losses and Maximum Pitting 


The data of table 3 indicate that the high-alloy 
cast iron E was considerably more resistant to cor- 
rosion, as measured by both weight loss and maxi- 
mum pitting, than either the plain cast iron A or the 
low-alloy cast irons B,C, D, NC, and N. It is note- 
worthy that, except in the deepest pit 
measured on the specimens of material E in more 
than 14 vears of exposure was only 74 mils. 

The effect of composition on the corrosion of the 
low-alloy cast irons can be observed to somewhat 
better advantage by calculating the weight losses of 


cinders, 
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these specimens on a relative basis. ‘This would have 
the effect of eliminating the differences in corrosivity 
of the soils due to differential aeration and emphasiz- 
ing the effect of graphitic corrosion. In table 4 the 
weight losses of materials B, C, and D are shown for 
each soil and for each period of exposure relative to 
Except for the well-drained 


the plain cast iron A 
and 62), and soil 


soils of high resistivity 
66. in which accelerated corrosion of the alloy cast 
irons occurred, the relative weight losses of the speci- 
mens during the initial period of exposure is seen to 
content of nickel 
vears’ 


soils 


eo e&Ff 
soils oh, eo 


Increasing 
the first 


losses of 


decreased with 
For example, during period of 5 

duration the weight material D in 
58 and 63 were only 20 and 16 percent, respectively, 
of the weight losses of plain cast iron in the same soils 
for the same period. However, this initial advantage 
from the addition of nickel was usually not main- 
tained, with the result that the weight losses of the 
alloy east irons exposed for the maximum period 
usually did not differ greatly from the losses of plain 
cast iron Hence, it would appear that the rates of 
corrosion of the alloys containing the higher amounts 
of nickel decrease less with time than do the rates for 
the lower amounts of nickel 


have 


the alloys containing 
and plain east iron 


5. Strength of Corroded Cast-Iron Pipe 


Corrosion of cast iron in soils is characterized by 
the conversion of the metal into a lave! consisting 
chiefly of oxides of iron and graphite. Although the 
original shape and appearance of the metal are re- 
tained, visual observation gives no indication of the 
It is generally recog- 
retains of its 

to which cast-iron 


extent of corrosion (fig. 1 

nized that corroded iron 
original strength, but the extent 
pipe may corrode underground and still retain suf- 
ficient strength to withstand the pressures commonly 
used in water- and gas-distribution svstems has not 


cast some 


previously been estimated 

In order to evaluate the residual strength of cor- 
roded cast-iron pipe, the samples of pipe that had 
been removed from the more soils after 
the longer periods of exposure were subjected to 
hydraulic The screwcap on one end of 
each sample was replaced by a similar cap in which 
a fitting had been inserted and the pipe section con- 
nected by copper tubing to a hand pump of suitable 
capacity. The pressure was increased at the rate 
of approximately 10 ('b/in.*)/see until failure of the 
pipe occurred the maximum pressure of 500 
lb/in.? was attained. After the application of hy- 
draulic pressure, the corrosion products were re- 
moved, and the condition of the specimens was 
evaluated in the usual manner. The hydraulic 
pressures applied to the samples and the number of 
perforations in the samples after being cleaned are 
given in table 5 

The data in table 5 show that most of the speci- 
mens withstood a maximum presssure of 500 Ib/in.’, 
although removal of the corrosion products revealed 
numerous holes of various diameters. The few low 


Corrosive 


pressures 


or 
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values of bursteng pressure given in the table 
ably should not be taken as accurately meas 
the strength of the corroded pipe because of « 
tion and probable damage to the pipe in har 
and shipping. Hence it is reasonable to 
that cast-iron pipe in an advanced stage of gra) 
corrosion is able to withstand the minimum pri 
required of class 150 pipe. 

In considering the data of table 5, it is impx 
to note that the pipe samples were subjected 
hydraulic pressure from 6 months to a vear 
removal from the test sites. Inspection of thi 
ples showed that during this interval some 
tion of the corrosion products occurred, as 
indicated by slight increase in the volume of 
corrosion products. Because oxidation is associat 
with loss in strength, it is reasonable to assume that 
the samples had not been removed from th 
sites, they would have withstood the applied 
sure, even if they had been in a more advanced 
of corrosion. Instances of the loss in strengi 
graphitized cast-iron objects on exposure to air ha 
been cited by Speller {4}. 

In order to evaluate the strength of the corrosi 
products of cast-iron pipe undisturbed by remoy 
from the test sites, 3-ft sections of 6-in. cast-in 
pipe, closed at both ends, were buried at 4 of 
Bureau test The investigation was 
taken with the cooperation of the Cast Iron P 
Research Association. Provision was made for a 
plying hydraulic pressure to the pipe sections 
means of a fitting inserted in the closure at one e1 
of the pipe According to the plan of the test 
pressure of 400 lb/in.? is applied to the pipe section: 
at approximately 2-vear intervals until rupture o 
curs, at which time the section will be removed a: 
the extent of corrosion evaluated. As none of 
pipe sections has failed vet withstand | 
applied pressure, the test is necessarily incomplet 
However, the data obtained to date in certain 
the provide significant information 
strength of corroded cast-iron pipe under 
conditions of exposure. 

In table 6 are shown the number of years duri 
which the pipe sections have withstood the appli 
pressure. For comparison, the condition of similar 
exposed specimens of cast-iron pipe, as indicated | 
the maximum and average penetration, is also show! 
The data reported for site 70 is especially significa 
Although 1 of the control specimens was perforat 
within 6 years, the 3-ft section has withstood t 
applied pressure for 24 years, during which time 
would be expected to be converted almost entir 
to corrosion products. Similarly, the developmer 
of a pit 302 mils in depth in 6 years on the cont! 
specimens in soil 118 indicates that the 3-ft sector 
which withstood the applied pressure after 24 year 
of exposure, should also be in an advanced stag 
of corrosion 


col 


Stag 


sites. unde 


as to 


on 
norn 


soils 


6. Summary 


This report contains the results of measureme! 
of weight loss, maximum depth of pitting, am 





c bursting pressure made on samples of probably to be considered minimum values The 
ist irons after exposure to different soil con- | results of incomplete field 


tests of the bursting 
Additions of nickel up to 3 percent signifi- | pressure of undisturbed sections of plain cast-iron 

reduced the initial corrosion of cast iron in | pipe are in good agreement with the results of the 

drained soils of low resistivity, but this measurements made on the stored samples 

age Was not maintained for the duration of 

osure tests. Samples of austenitic cast iron 

is 15 percent of nickel together with 


The direction of this program by Mr kK H 
Logan prior to his retirement in 1946 is acknowl 
edged 


m and copper, showed good resistance to 
yn in most of the soils, especially with respect 
mum depth of pitting 
some exceptions, samples of severely cor- 7. References 
lain and low-alloy cast-iron pipe withstood a 
_ W. A. Wesley, H opece 
im hydraulic pressure of 500 Ib/in? after & Allova 2° 325 (19: 
e at the test sites for periods up to 11 years 2] Irving A. Denison and M« 
rage in the laboratory for approximately 1 NBS 44, 47 (1950) RP2057 
Because the effect of oxidation in the air is : wy _ ae ometen, TS Shp 
7 OU (1945 
Ken the laver of corrosion products, the | F. N. Speller, Corrosion causes and preventio a 
d values for hydraulic bursting pressure are Hill Book Co.. New York. N. ¥.. 1951 


vin Romanoff 


TABLE 1 Properties of the soils at the test 


Percent Ohm-cm 
Atlanta, Ga ; 3.7 ) 17, 800 
Loch Raven, Md 32.0 21 
Meridian, Miss 4.6 6. 920 
Wilmington, Calif 2.4 
Phoenix, Ariz 16. 5 7% 232 
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‘alorimetric Properties of 41° and 122° F Polybutadienes 
George T. Furukawa and Robert E. McCoskey 


f 41 
e Was to 


he 122° } polymer 


to be 200° to 295° K an 


ithed tab 


rhe data were ‘ to construct a smo« 


from O° to 330° hk 


1. Introduction 


metric investigations of the thermal proper- 
1° and 122° F butadiene-styrene copolymers 
ning 8.58 weight percent of bound styrene are 
ed in a previous paper [1 The properties 
d included heat capacity, enthalpy, entropy, 
transformation, and crystallization The two 
which were prepared at two different 
ratures, exhibited wide differences in thermal 
ties. For example, the glass-transformation 
erature of the 41° F copolymer was found to be 
K and of the 122° F copolymer, 193° K. Also, 
F copolymer exhibited crystallization from 
to 285° K, whereas the 122° F copolymer 
d no such crystallization effects. This paper 
vith similar calorimetric studies with 41° and 
Fk polybutadienes, both of which have been 
to undergo ery stallization The effect of the 
of crystallinity upon the heat capacity of the 
volymers was investigated and the degrees of 
llizability of the two polymers were compared 
heat-capacity Measurements made from 
16° to 330° K. The results were used to com- 
a smoothed table of heat capacity, enthalpy 
ntropy from 0° to 330° K 


2. Apparatus and Method 


Viners 


were 


details of the calorimetric 
«l used have been described [1, 2 
was sealed in a copper container within an 
atic calorimetric system. The container was 
led with a platinum resistance thermometer, a 
and a vane system for rapid distribution of 
To thermally isolate the container during 
lorimetric experiments, the temperature of the 
system surrounding the container was main- 
the same as that of the container surface by 
s of differential thermocouples and shield heat- 
The surrounding the container was 
iated, and the container and the adjacent shield 
es were polished. The electric power input 
neasured by using a Wenner potentiometer, and 
ngth of the heating periods was measured by 


of a timer operated on standard 460 eveles 
— ? 
rk discussed herein was 


apparatus and 
The polymer 


space 


performed as a par f the research project 
the Reconstruction Finance Corporat { Synthetic Rub 
tion with the Government synthetic rubber program “a7 


n brackets indicate the literature references at the end of this paper 


uutadienes was meas 


ired from 15° to 330° K 
er (195° K) glass-transformatior 
} 


polymer exhibited a gher 


ten pera 
legree 
allization temperature range of the 
ymer 200 270° Kk 
thalpy, and entropy 


The temperature measurements above 90° K are in 
accordance with the 1948 International Temperature 
Scale [3]; between 10° and 90° K they are On a provi- 
sional scale [4], which consists of a set of platinum 
resistance thermometers calibrated with a helium-gas 
thermometer. 

The polymers were subjected to two generally 
different rates of cooling prior to the heat-capacity 
experiments in order to study in what way their heat 
eapacity and ery stallinity were dependent upon the 
thermal history. In one series of experiments, the 
polymers were cooled as rapidly as possible by im- 
mersing the calorimeter in liquid nitrogen with 
helium gas in the space surrounding the container 
By this procedure the polymers were cooled to about 
90° K in 30 min and to 80° K in an additional 30 
min, and a portion of the polymer molecules was ex- 
pected to be “frozen” in the higher energy states 
In the second series of measyrements, the polymers 
were cooled slowly in order to obtain as high a degree 
of crystallinity as practicable. By maintaining a 
high vacuum in the space surrounding the container 
and by using different refrigerants successively, the 
cooling process was prolonged from a few days to 12 
days. 


3. Samples 


The two polybutadiene samples were obtained 
through the courtesy of the University of Akron 
Government Laboratories, Synthetic Rubber Re 
search. The 41° F polybutadiene designated MS 
1045 was prepared using a low-sugar cumene hydro- 
peroxide-redox recipe [5], and the 122° F polymer 
designated GL-657 was prepared by using a recipe 
[6] generally used for GR-S type polymers. The 
polymerization formulas are given in table 1. 

The polymers were purified by M. Tryon, of the 
Rubber Section of the Bureau, by a_ procedure 
previously outlined [1]. The analysis of the samples 
is given in table 2. The samples were pressed into 
sheets \ to \ in thick and cut to fit between the 
vanes of the calorimeter containers Before sealing 
the containers, the polymers were pumped at high 
vacuum for 3 days to remove moisture, air, and 
benzene. The mass of the MS-1045 sample in- 
vestigated was 42.706 g, and that of the GL—-657 
was 47.081 g. A small quantity of helium gas was 
sealed with the polymers to enhance thermal equi- 
librium during heat-capacity experiments. 
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Heat treatments and the observations with th 
polybutadie ne, VS 1045 


Temperature of 
J ; Temp. 

drift ot 
y 


range 
r r f " » ake O 
eatment of samples measure- 


ments 


lly from room temper 


room temperature for 1 day 
from room temper 
K over 2 days 


from room tem 
and then to 15° K 

from room temper 
kK 


per 


and then to 53 


Observed heat capaciti« 
VM S-1045 


4. Results 
4.1. 41°F Polybutadiene, MS-1045 


The experimental details of heat treatment, 
temperature range of the measurements, and tem- 
perature drifts with the 41° F polybutadiene are 
summarized in table 3. The observed heat capaci- 
ties are given in table 4, and the general behavior of 
these results are shown in the plot of figure 1. The 
measurements of run 1 were made after shock-cooling 
the polymer from room temperature to 77° K experiments, these molecules acquired _ sufficier 
During the course of the heat-capacity measure- | thermal energy for transition to lower energy states 
ments, upward temperature drifts were observed | with the consequent liberation of heat. Prior 
from 177° to 283° K as the polymer molecules slowly | run 3, the polymer was cooled slowly from roo! 
transformed into states of lower energy [1]. During | temperature to 77° K over 2 days, and the results 
the shock-cooling process a certain portion of the | showed upward temperature drifts from 18% 
polymer molecules failed to undergo transition to | 280° K. Below the glass-transformation temper 
states of lower energy. When the temperature was | ture (195° K), the heat-capacity results of run 
raised to a certain range during the heat-capacity | (polymer slowly cooled) were slightly lower than ! 
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if capac 


WS 


l This indicates more melting than crystallization ; 
consequently, the upward drifts observed are prob- 


ably the result of the formation of new crystals and 
of the recrystallization of polymer near the heater 
where the temperature was higher during the heating 
period. The broken line of figure 1 is the heat capac 
itv of hypothetical amorphous polymer and was 
obtained by linearly extrapolating the heat capacity 
in the range 295° to 330° K to the glass-transforma- 
tion temperature. This linear extrapolation of the 
amorphous portion of the heat capacity was made 
on consideration that the heat capacity of HyCar 
O. R.-15 [7], which does not crystallize, could be 
represented within 0.2 percent from 250° (just above 
the glass-transformation temperature) to 340° K by 
a linear equation. Similarly, in another noncrystal- 
lizing synthetic rubber, GR-S [8], the heat capacity 
from about 225° to 330° K could be represented by a 
quadratic equation to 0.1 percent or by a linear 
equation to about 0.3 percent or better 
The data of table 4 were used to construct the 
heat-capacity values of table 5, spaced at equal 
temperature intervals. The heat capacities given 
were based on the results with the polymer slowly 
cooled From 180° to 305 K the heat-capacity 
values of table 5 were obtained from al large scale 
plot of the observed heat capacities In the erys 
tallization interval (200° to 295° K) the results 
from the broken line were tabulated. Below 180 
K and above 305° K the deviations of the observed 
heat capacities from empirical equations were first 
| lotted, then the deviation curves and the empirical 
equations were used to obtain heat capacities at 
even temperatures The heat capacities below the 
lymer rapidly cooled), with the difference | lower limit (15° K) of the measurements were 
ising with temperature to a maximum of about obtained by fitting a Debye function to the observed 
ent near the glass-transformation temperature values in the temperature range 15° to 30° K. The 
ipward temperature drifts from about 200° to | function used is 
K are attributed to slow crystallization of the 
ner molecules. Except for a few heat-capacity 
s of run 1, all experimental points in the crystal- 
yn interval were above the broken line of figure 


323 


0.2425 D(88.7/T 





The glass-transformation temperature of this par- | Taste 6. Heat tr uments and th observation 

ticular sample of 41° F polybutadiene was found to polybutadiene, GL 

be about 195 iN Above 295 K the polymer Is 

believed to be completely amorphous 
The enthalpy and entropy of the 41 

obtained by evaluating the thermody- 


F polybuta- 


diene were 


imic relations 


I] 


respectively where H Hox and S So°K 
ire the enthalpy and entropy of the polymer relative 
to the absolute zero of temperature The other 
symbols have their usual significance. Except be- 
low 15° K and between 180° and 305° K eq (2) and 
;) were evaluated by tabular integration, using 
four-point Lagrangian integration coefficients [9] TABLE 7 
Below 15° K the Debve function (eq 1) was evalu- 
ated analvtically Between 180 and 305 K the 
enthalpy was obtained directly from the experimental 
measurements [1] by summing the energy input of 
run 3 (polymer was cooled slowly), in which the 
experiments were made continuously from 180.0093 
to 3OLGOISS kK When corrected to the even tem- 
perature interval (180° to 305° K), the enthalpy 
change amounted to 247.84 abs j g The entropy 
change was obtained by summing similarly the 
various (A7/7,,’s of run 3, where CAT is the en- 
thalpvy change of the heating interval and 7, is the 
midtemperature of the interval After correcting 
to the temperature interval 180° to 305° K, the 


oO 


entropy change amounted to 1.018 abs j °K 


4.2. 122° F Polybutadiene, GL-657 


The experimental results with the 122° F poly- 
butadiene are summarized in tables 6 and 7. Table 
6 shows details of the heat treatments, temperature 
range of the measurements, and the drift observa- 
tions. The observed heat-capacity values are given 
in table 7 and plotted in figure 2 When the poly- 
mer was cooled rapidly, the subseq uc nt heat-capacity 
experiments (run 3) showed upward temperature 
drifts from 148 to 253 kK However when the 
polymer was cooled slowly, the heat-« apacits results 
run 2) showed downward temperature drifts from 
186° to 189° K and upward drifts from 220° to 
265° K. These results show that in the experiments 
with the polymer rapidly cooled, the upward tem- 
perature drifts set in about 40 deg below the tem- 
perature at which the downward drifts began with 
the polymer slowly cooled. In the experiments 
with the 41° F polybutadiene slowly cooled over 2 
days, no downward temperature drifts were observed 
below the glass-transformation temperature. Only 
upward temperature drifts were observed. Prob- The upward temperature drifts from 200° to 270 
ably 2 days were not sufficient for much of the | K are considered to arise from slow crystallization o! 
41° F polymer to transform to lower energy states | the rubber polymer in a manner similar to the 4! 
at temperatures below the glass-transformation | F polymer described earlier. Also, the relat vely 
temperature high heat capacity above the broken line indicates 
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I22°F Polybutodiene 


capacity (see fig. 2 A large amount of crystalliza- 
tion has caused the apparent heat capacity to be 
relatively low 

The results of these measurements were used to 
compute a table of heat capacities at equally spaced 
temperature intervals (see table 8 The results 
used were those from the experiments with the 
polymer slowly cooled, and the procedure used in 
obtaining the smoothed heat capacities, including 
those in the crystallization interval, was similar to 
that with the 41° F polymer. The Debye function 
for the extrapolation of the data below 15° K was 


0.2631 D(89.5 


The glass-transformation temperature of this ma- 
terial was found to be 187° K. The erystallization 
range was 200° to 270° K, which is about 25 deg 
narrower than the 41° F polybutadiene (200° to 
295° K Above 270° K, where the heat capacity is 
smooth, the polymer is believed to be amorphous 
The enthalpy and entropy of the 122° F poly- 
butadiene were evaluated by a procedure similar to 
the 41° F polymer. The values are given in table 8 


4.3. Reliability of the Results 


The reliability of the heat-capacity measurements 
with these polymers is difficult to evaluate because of 
the nonreproducibility of the physical state from 
experiment to experiment. The heat-capacity meas- 
urements made previously with normal materials of 
low molecular weight are believed to be accurate to 
0.2 percent [10]. The results with these polymers 
show considerable scattering i the temperature 

melting than crystallization. In the experi- | range where temperature drifts occur, as the results 
with the polymer slowly cooled (run 2) all the | are dependent upon how long temperature equilibri- 
s were above the broken line. However, with | um is awaited. In the experiments where tempera- 
lymer rapidly cooled (run 3) a low point was | ture drifts occurred, although no attempt was made to 
ed, which gave a y-like shape to the heat | wait until complete temperature equilibrium was 
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» on . 
| TABLE 0 ( om parison of the enthalpy changes of 41° and 
| : 


reached. the general nature of the heat capacity ts 
believed to be preserved In the temperature range 
where no temperature drift occurred, the heat- 
capacity results obtained for a given heat treatment 
are believed to be accurate within 0.2 to 0.3 percent 

The enthalpy and entropy values given in tables 5 


and S& are believed to be generally accurate to OQ.o 


polybutadienes from 175° to 305 


pe reent 


5. Discussion 


In the investigation of the 41° F butadiene-st 


The results of these experiments show that by i : 
copoly mer [1] of S.08 weight percent of bound st 


changing the polymerization temperature of polvy- : 
ares, ed - - to 199 KF aoa ass slic cn. the enthalpy change difference of the anneales 
tion temperature was lowered from 195° to 187° K integrated results amounted to about 16 to 19 al 
This behavior is similar to that observed with 41 J & In this Investigation the 41° and 12 
and 122 F butadiene-stvrene copoly mers containing poly butadienes gave 33 and | 4 abs j £ respect ‘ 
$.58 weight percent of bound-styrene [1] previously Bekkedahl and Matheson [11] reported 16.71 
reported. Also, the temperature range of erystal- | ® for the heat of fusion of natural rubber, wi 
lization has become narrower in the 122° F polymer indicates that 41° F polybutadiene has a hig 
200° to 270° K The crystallization temperature erystallizability than natural rubber. 
range of 41° F polymer was found to be from 200° to 
295° Kk In the butadiene-styrene copolymers 1), The authors are indebted to W Bruenner, . 
crystallization took place in the 41 F copolymer but | King. and M. L. Reilly for assisting with som 
was absent in the 122° F copolymer nt the experimental work and computation, to M 
— _ ° th age changes irom ate wo Tryon for purifying the samples, and to R. A. Pa 
105° K are given for both polymers The results | con for the analvsis. 
listed are from experiments in which the polymers 
were given various prior heat treatments, as described 6. Referen 

. ces 


earlier. Slight adjustments have been made in the 


original experimental enthalpy in order that the [ rr. Furukawa, R. E. MeCoskey, and G. J. King 
Research N BS $0, 357 (1953) RP2425 
yy 3. Scott, C. H. Meyers, R. D. Rands, Jr., 1 
n the 41° F polymer the enthalpy changes betwee , _— . 
In the amy: tae a ‘. _—— wo Brickwedde, and N. Bekkedahl, J. Research NBS 35, 
shock cooled and annealed differ by only 5.51 abs 39 (1945) RP1661 
jz In the case of the 122° F polybutadiene, the [3] H. F. Stimson, J. Research NES 42, 209 (1949 
difference is 12.07 abs JZ The enthalpy changes 4) H J Hoge and F. G. Brickwedde, J. Research 
obtained for the 41° F polybutadiene for this temper- mW pu RI 1188 a ald a 
. ate communicatioir rom ‘ t e 

ature interval are much greater than that of the 122 a Neng; Sear Lahacttesion to L 
I polybutadiene These results seem to indicate National Bureau of Standards 
that the 41° F polymer has a higher crystallizability R. W. Laundrie, Butadiene-styrene copolymer 
than the 122° F polymer. The integrated enthalp refractive index studies, Project P 100.74, AU-725 

| egrated en hal , l N. Bekkedahl and R. B. Scott, J. Research NBS 29, 87 
change was obtained by tabular integration of the 1942) RP1487 
heat capacity given by the hypothetical amorphous | R. D. Rands, Jr., W. J. Ferguson, and J. L. Prather 
curve. The comparison of the enthalpy change in _ Research NBS 33, 63 (1944) RP1595 
the 122° F polymer when sl k-cooled wit . [9] Tables of Lagrangian Interpolation Coefficients (( 
ntevrated al ’ i ‘at in tl the pone | os & bia University Press, New York, N. Y., 1944 
integrated Value ine ea es mt the amorphous con- 110) G. T. Furukawa. D. C. Ginnings, R. E. MeCoskey. a 
dition was practically frozen-in On the other hand, R. A. Nelson, J. Research NBS 46, 195 (1951) RP2 
a similar comparison of the 41° F polymer would [ll] N. Bekkedahl and H. Matheson, J. Research NBS 1i, 

503 (1935) RP844 


temperature range be the same in both polymers. 


indicate that even upon shock-cooling considerable 
crystallization took place | WasHINGTON, July 30, 1953 








Thermal Degradation of Tetrafluoroethylene and 
Hydrofluoroethylene Polymers in a Vacuum 


S. L. Madorsky, V. E. Hart, S. Straus, and V. A. Sedlak * 


refion and tetrafluoroethvlene photopolymers, on pyrolysis in a vacuum at 423.5 
513.0° C, vield almost 100 percent of monomer The rate of formation of monomer at a 
given temperature follows a first-order reaction and is independent of the method of prepara 
tion of polymer or its initial average molecular weight. The activation energy was deter 
mined by a pressure method and a weight method, and a value of 80.5 keal was found by 
both methods reliminary heating of Teflon in air at 400° to 470° C did not chang 
appreciably its rate of degradation into monomer when it was subsequently heated in a 
vacuun Polvviny ioride, 1,1-polyvinylidene fluoride, and polytrifluoroethylene wer 
pvrolyzed in the range 372° to 500° C rhe volatiles consisted in all cases of HF and a 
wax-like material consisting of chain fragments of low volatility Polyvinyl fluoride and 
polytrifluoroethvlene degrade to complete volatilization, whereas 1,1-polyvinylidene fluoride 
becomes stabilized at about 70-percent loss of weight he rate-of-volatilization curves indi 
cate a first-order reaction for polyvinyl fluoride, a zero-order reaction for trifluoroethylene 
and an undetermined order for 1,1-polyvinylidene fluoride. The order of thermal stability 
for these polymers, as compared with polymethylene, is as follows: Polyvinyl] fluoride 
polyvmethvlene polytrifl oroethvlene ] l-polv vinvlidene fluoride polytetrafluoroet hvlen 


l er I 8 hyudrofiuorocarbon poliyme 


1. Introduction Papi 


is very little in the literature on the thermal] 
ation of fluorocarbons. in veneral, or on 
irbon polymers, in particular. Swarts [1],’ 
sand Cady [2], and Steunenberg and Cady [3] 
zed a number of low molecular weight fluoro- 
is in the presence of a glowing platinum 
it. Lewis and Naylor [4] pyrolyzed polytet- 
ethylene at 600° and 700° C and at pressures 
ng from 5 to 760 mm Hg The volatiles consisted 
Fy, C3F,, and C,F,. In this work a study was 
of the thermal degradation of a series of 
polymers to determine their relative thermal 3A . 
aratus and Experimental Procedure 
ty, the nature and relative amounts of the PP P 
es given off, and the rates of thermal degrada- The investigation of thermal degradation of this 
This series includes polytetrafluoroethylene | series of polymers was carried out along two lines 
n C, F, ® polyy invl fluoride | Cc H F ® 1. Pyrolysis in a vacuum and fractionation and analysis of the 
lvvinylidene fluoride [—-C,H,F and poly- | velatile products. This procedure was followed for all of the 
oethvler : C.HF a ' polymers, except Teflon, using a Dewar-like molecular still, 
: uM : " which has been described in previous papers 6, 7] A 20- to 
{0-mg sample, either in solution or in finely divided from, was 
: spread on a platinum tray. Size of the sample was limited so 
2. Materials Used as to prevent loss of material by spattering during pyrolysis 
The sample was first subjected to a preliminary heating in a 
polytetrafluoroethylene was a commercial vacuum for 2 hr at about 150° C in order to eliminate the 
tape, 0.07 mm thick The polyvinyl fluoride solvent and adsorbed gases It was then brought to the 


: ; , 3 : ; ; ‘ temperature of pyrolysis by heating for 45 min, and this 
ie polyvinylidene fluoride were prepared by | | mperature was then maintained for 20 min. The following 


du Pont de Nemours & Co lhe polytri- fractions were collected: I, residue; II, a waxlike material 
f thy lene was prepared from the monomer by nonvolatile in a vacuum at room temperature; III, a fraction 
olymerization at 20° C in the presence of | volatile at room temperature _IV, @ gaseous fraction noncon- 
butyl peroxide and then heated overnight at Soaines ae Oe Senegal < e~te-v + ; ae = 
. ’ . four fractions were determined, and, in the case of fractions 
‘ Analyses for ( ? H, and k in the hydro- II] and I\ . chemical com positior Was ce t¢ rmined by means of 
arbon polymers are given in table 1 the mass spectrometer To facilitate mass-spectrometric 
——— , analysis, fraction II] was further subdivided by distillation 
iress: U. S. Public Health. Atlanta. Ga into a light fraction, IIIA, and a heavy fraction IIIB. For 
k Was performed as a part of the research pr t hig temperature Teflon, which requires a higher temperature of pyrolysis and 
: — Devree wey be ~ a gdh ig bam ~ wader Art which yields almost 100 percent of monomer, a different twpe 
nistry Division, September 19% of apparatus was used 
brackets indicate the literature r ence t the end of this pay 2. Rate of volatilization of polymers na vacuun This 
mer and polymer were prepared by R . = Be . property was investigated in the case of all polymers by a 


mer Structure Section of the Natior = 
ps weight method, and in the case of Teflon, also by a pressure 


prepared by the method of Park 





method The spring-balance appa d the weight Rates of volatilization of Teflon were determing 


method is described in a previous paper [8 The sample was the pressure developed in the calibrated volume 
d to 5 to 8 mg in order to avoid loss by spattering during pressure Was measured at time intervals by means of a 


g. It was placed in : ible, and the plying manometer, F, containing a low-vapor-pressure s 
le was suspended by means of a 3 tungsten wire from oil in the left arm of the manometer on top of a n 
ingsten-spring b ! ‘ 1 in a Pyrex column The manometer was calibrated by means 
ich could be evacuat » about | mm Hg three-scale McLeod gage, reading to 25 mm Hg, 
icible was heated externally, and loss « ght of the precision of 0.2 percent The scale on the mult 
was determi! ' cre e on an exten manometer could be read with a precision of 0.02 
Pressures up to 6 mm were meas ired 
apparatus fo tud ‘ of thermal degradation : A sample of the fraction noncondensable at the temper 
Rone ammatically of liquid nitrogen, corresponding to fraction IV whe 
} em long, Dewar-like apparatus is used for pyrolysis, was obta 
closed | the following manner. Liquid nitrogen was placed a 
tube, B, trap L The condensable material, corresponding to fr 
ind joint, | IIT, condensed in this trap, while fraction IV remaine 
from 5 to | pended in the space between pump E and stope 
ng occurre d Stopeock H over t ibe G was closed and the tube was 
aterial in off at K The contents of tube G were analyzed in th 
the study spectrometer The weight of fraction IV was deter 
from its pressure, total volume, and analysis 
ire. Was lo obtain a sample of fraction IIT for analysis, the s 
center for | was evacuated while trap L was still immersed in liquid 
noved in gen With stopeock D closed, fraction IIT was then a 
approxi to expand into the apparatus by removing liquid ni 
ec 110-v | from around trap L. Fraction III was then sampled i 
stabilizer | G’ by closing stopcock H’ and placing liquid nitrogen a 
perature G’ while the sample tube was sealed off at K We 
inum-1l0 fraction III was determined from its pressure, total 
i i and analysis 
moved in In some experiments it was found expedient to coll 
nocoupk entire fraction III and to weigh it his was done by 
tem ing it in one of the weighed tubes, M, provided with 


perature inside of tu A was calibrated ans of ar joints, O. Liquid nitrogen was placed around the lower 
additional platinum versus platinu reent lium of M and held there until condensation was complet« 
ermocouple placed temporarily) le of 1 tu here tube was then sealed at N, without melting it off. and 
appreciable temperatu gracdit t tl l l the sequently weighed tepeated experiments showed a 
Ten perature iat I vot hermo- check betweer he two methods of determini! g total 
about 0 deg of fraction ITI 
, Was The residue, fraction I, was weighed in tube A A wa 
apparat. vacuat o about material, fraction Il, appeared in some experiments as 
ns of an oil imp. not ) in the posit in the cold part of tube B, just outside the heater 
igure, and a | itrogen trap r ating of a sample weight of this fraction was determined by subtracti 
of Teflon in a vacuum at about 300° t for 48 hr | sum of weights of fractions III and IV from the tota 
sulted in no appreciable loss in \ f he ‘ _ which weight of the sample In cases where fraction II d 
been maintained overnight at t oO , ten ‘ . appear, there was a good balance between total weight 
then moved into position aro t t took 5 fractions and original weight of sample 


j for the tem perat ure of t< ach a constant 


value At the termination of ar ! B. f ater Was 4. Pyrolysis of Teflon 


removed quickly from the apparatus ; 

During pyrolysis the apparatus was cut o the | Results of py rolysis of Teflon are shown in table 2 
evacuation pump by means of stopcock while a met - | This table gives also, in the last column, the result 
Tusi pump, FE, 1 —_ pera of studies of rates of thermal degradation. The « 
wzainst a back pressure of 25 mn periments were carried out in the pressure apparat 


effective in removing the volatiles fror 
this pump and stopcock D was calibrate shown in figure | and are marked P in column 


vec 


¢ study of rates of thermal dé qg adation of polyme S 


hy the pressure method 
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fraction I] deposit 
Fraction iV, on the 
amounted to 0.1 mole percent of the total 
ed part and consisted of CO Thus fraction 
ractically equal to total loss in weight viven 
table 2 
about 50 


The re Was ho wax-like 
of these experiments 


mn 5 of Teflon retains its original 
intil volatilization 
this point it softens and slumps 
thod of comparing the thermal stability of a 
of polymers is described in previous papers 
In this work the thermal stability of Teflon 
n in figure 2, compared with that of hydro- 
irbon polymers and polymethyvlene. In the 
leflon, the time required to heat a sample to 
temperature only 3 ta 5 min 
pared with 45 min for the other polymers 
re pyrolyzed in the Dewar-like apparatus 
to put the thermal-stability curve for Teflon 
ymparable basis with the other curves, it was 
15° C to the left This adjustment was made 
experiments with polystyrene and 
fluoroethylene in the Dewar-like apparatus 
pressure apparatus 
ne experiment a sample of about 100 mg of 
was heated in four consecutive steps to almost 
te volatilization. In large size 
sample, the volatiles were condensed during 
sis In trap Q, figure 1, by means of liquid 
n to insure their complete removal from the 
one. However, when collecting samples for 
for analysis, or evacuating the 
itus between steps, trap Q was maintained at 
temperature. After each step, samples from 
ms Ill and IV were collected for analysis in 
G and G’, respectively fig. 1 
Ill was obtained by collecting it in one of 
bes M and weighing. Allowance 
nount of fraction III collected in tube G 
s of the step-experiment are shown in table 3 
on II appears only in the last step, where the 
itive volatilization was 93.7 percent. In all 
essure experiments shown in table 2, except 
t one, maximum volatilization was 76 percent; 
iction IT was not observed in these experiments 


percent ot 


erating was 


basis of 


view ol the 


ng or when 


Total weight of 


was made for 








In the last experiment, at 513.0° C, where volatiliza- 
tion was 90.1 percent, the weight of the sample was 
small; and, if a fraction II appeared, it was too small 
to be detec ted 

Results of mass-spectrometer analysis are 
shown in table 3. The CO shown in this table ap- 
peared as fraction LV, and the rest appeared as frac- 
tion III. As seen from this table, composition of the 
volatiles did not vary up to at least 71.4 percent of 
volatilization. The Sify, CO., and CO might have 
resulted from oxygen in the polymer or, most likely, 
from a reaction between Teflon and quartz [10 


also 


5. Rates of Volatilization of Teflon 


5.1. Determination by the weight method 


data on rate with 
the spring-balanc« 
shown in table 2. The weight experiments are 
marked W in the first column of this table. In 
ficure 3 the rates are plotted in percentages of the 
original sample per minute against percentages vol- 
atilized. By extrapolating the straight lines to zero 
volatilized, the initial obtained. These 
initial rates are shown in the last column of table 2 
for experiments marked W 


experiments 
apparatus are 


Experimental 
Teflon tape in 


rates are 
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of thermal degradation of 
entage 


Fieure 3. Rate 


method, as a function of pe 


5.2. Determination by the pressure method 


Inasmuch as Teflon yields, on pyrolysis, mainly 
the monomer and a small percentage of other prod- 
ucts, all in constant proportion and volatile at room 
temperature, the rate of volatilization can be ob- 
tained from a curve in which pressure of the gas in a 
fixed volume is plotted against time. Results of a 
series of rate experiments in the pressure apparatus 
are shown in the last column of table 2. An example 
of a pressure-time curve is shown by the curved line 
in figure 4 for a rate experiment made at 513° C 
This experiment is shown in table 2 as the last 
experiment of the series 

The condition for a first-order reaction is that the 
plot of In(a—z), where (a—z) is the residue, against 
time, ¢, is a straight line [11]. In such a plot, rate of 
reaction, &, is the slope of the straight line, and can 
thus be evaluated. The straight line in figure 4 was 
obtained by plotting log of residue, calculated in 
terms of pressure, against time for the 513.0° experi- 
ment. Using such curves as shown in figure 4, the 
rates shown in the last column of table 2 were ob- 
tained for the experiments marked P 

To determine the effect of heating Teflon in air on 
rates of volatilization, samples were heated in air at 
400°, 425°, and 450° C. The amounts volatilized 
were about the same as in a vacuum. When the 
samples were heated at 470°, the amount volatilized 
was slightly greater in air than in a vacuum. The 
residues from air-heated samples were pyrolyzed in a 
vacuum. The behavior of these with 
regard to rate and products of volatilization, was 
about the same as for original Teflon 

The effect of thickness of the Teflon specimen on 
rate of volatilization in a vacuum was determined by 


residues, 











Rate of thermal degradation of Teflon at 51 
by the pressure method 


Ficure 4 


, Time versus pressure plot time versus log of residue | 


the following experiment. A’ 1.3-mm_ cube of 
Teflon was pyrolyzed at 514.5° for 133 min. Tota 
volatilization was 78.3 percent, and the rate was th 
same as for tape Teflon pyrolyzed under simila 
conditions. In another experiment, a 1.9-mm culy 
of Teflon was pyrolyzed at 495°. Initially, the rat 
was slow, but it soon reached the value for tap 
Teflon. 

Some photopolymers of tetrafluoroethylene wer 
prepared [12] by irradiating the monomer wit! 
ultraviolet light in the presence of (CF;),Hg, benzo, 
peroxide, CF;1, and di-tert-butyl peroxide as cata- 
lysts; also in the absence of a catalyst. On pyrolysis 
in a vacuum, these polymers had the same rates o 
volatilization as the tape Teflon. 

The activation energy is obtained by plotting 
logarithms of rate versus inverse of absolute tempera- 
ture. This plot is shown in figure 5 for the weight 
and the pressure experiments. All the points fall o1 
a straight line. Multiplying the slope of this line by 
2.303R, where R is the gas constant, 1.987 cal/deg 
value of 80.5 keal is obtained, which is the activatior 
energy for the thermal degradation of Teflon. Th 
frequency factor, A, as calculated by means of 
Arrhenius’ equation, has a value of 4.710", wher 
rates are expressed in fraction per second 


6. Pyrolysis of Polyvinyl Fluoride, Polyviny- 
lidene Fluoride, and Polytrifluoroethylene 


Pyrolysis of these polymers was carried out in th 
Dewar-like molecular still. Experimental details 
are shown in table 4. The gaseous fraction, I\ 
amounted in all cases to less than 0.1 percent of the 
total volatilized part. Mass-spectrometer analysis 
showed it to consist of hydrogen and carbon monox- 
ide. The less volatile fraction, I11, was found o1 
mass-spectrometer analysis to consist, in the cases 0! 
polyvinyl and polyvinylidene fluorides, mainly 0! 
SiF,, H,O, and some unidentified hydrofluorocarbor 
fragments. In the case of polytrifluoroethy!en 
this fraction, in addition to SiF, and H,O, contained 
some CQ). 
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Sif, and H.O are due to a reaction between 
esulting from pyrolysis, and SiO, in the glass 
atus. This reaction can take place at room 
rature in the presence of a trace of H,O For- 
n of H,SiF, from HF and Sif, is possible, but, 


331 


under conditions of vacuum and insufficient amount 
of H,O, hydrofluorosilicic acid tends to decompose 
back into HF and SiFy. Pvrex also contains some 
B.O;, which might react with HF, but there was no 
evidence in the mass-spectrometer analysis of the 
presence of BF Aluminum oxide, ordinarily pres 
ent in the glass, does not react appreciably with HF 
at room temperature 

The CO, found in fraction II] from polytrifluoro 
ethvlene could not be explained on the basis ot a 
reaction between fluorocarbon and SiO, of the glass 
According to White and Rice 10] such a reaction 
takes place at elevated temperatures, and in this case 
the polymer was in contact during pyrolysis with 
platinum, and the polymer fragments were in contact 
with glass at low temperatures. It is possible that 
the polytrifluoroethvlene contained some OXVgen as 
a part of its structure. This oxvgen, on reacting 
with the polymer, could give CO,. No analysis for 
oxygen in the original polymer was made because of 
the difficulty involved in such analysis in the presences 
of fluorine. The fluorocarbon tends to react with 
the glass container at the elevated temperature 
emploved in analysis, and CO, results from such a 
reaction. The CO detected in fraction IV from 
p\ roly SIS of hvdrofluorocarbon polymers could come 
from oxygen as an impurity 

Fraction III was calculated in the case of all three 
polymers to HF on the basis of Sif, content in the 
volatiles. This fraction is shown in the last two 
columns of table 4 in weight percent of the total 
volatilized part and in weight percent of sample 
Fraction I] consisted of a nonvolatile light-brown 
wax-like deposit, soluble in acetone. Because frac- 
tion IV was only about 0.1 percent, fraction I], in 
percentage of total volatilized part, can be taken as 
the difference between 100 and the percentage of 
fraction III, in terms of volatilized part, given in 
table 4. Mass-spectrometer analysis of fraction IT] 
indicated the presence of hvdrofluorocarbon mole- 
cules of molecular weight up to 150. However, in 
calculations of HF as fraction II], these hydrofluoro- 
carbon fragments were included in fraction II be- 
cause their nature could not be identified. Judging 
from results of pyrolysis of other polymers, in the 
temperature range of 370° to 430° C [7, 9], fraction 
I] from the hydrofluoropolymers should have an 
average molecular weight of about 600 to 700 The 
residue, fraction 1, from all three polymers appeared 
light brown during the early stages of degradation 
and dark brown toward the end. In the case of 
polyvinylidene fluoride, the residue, above 50 per- 
cent of volatilization, appeared black 

As can be seen from table 4, polyvinyl fluoride and 
polvtrifluoro ethylene volatilize almost 100 percent 
Polyvinvlidene fluoride, however, seems to become 
stabilized at around 65 percent of volatilization. 
When the temperature of pyrolysis was raised from 
156° to 530° C, the additional loss in weight was only 
5.5 percent. In one experiment in the pressure 
apparatus a 69.3-mg sample of polyvinylidene fluo- 
ride was heated in a vacuum from room temperature 
to 650° Cin 80 min. Total loss by volatilization was 
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The residue was further heated in the 
apparatus for 30 min by means of a Meker burner, 
applied to the outer tube B (fig. 1 Additional loss 
in weight was 1.5 percent of the original sample 
The residue resembled coke in hardness and appear- 
ance In another experiment a 78.9-mg sample of 
polyvinvlidene fluoride was heated gradually in the 
same apparatus trom temperature to 500 
during 2 hr and then kept at 500° for 1 hr. Loss of 
weight by volatilization was 70.3 percent. Chemical 
analysis showed that the residue still contained 1.7 
percent of H and 12.6 percent of F by weight 

Relative thermal stability curves for the 
hvdrofluorocarbon polymers are shown in figure 
in comparison with similar curves for Teflon and 
polymethvlene. The polymethylene curve is based 
on pyrolysis experiments made in the Dewar-like 
apparatus at five temperatures. Experimental con- 
ditions were about the same as in the case of the 
hydrofluorocarbon polymers. The polymethylene 
was a nonbranched high-molecular weight polymer 
of the same stock that was used by Mandelkern and 
13] in their study of intrinsic viscosity 
HF from the three hvdrofluorocarbon 
polymers during pyrolysis, in percentage of available 
HF in the polymer, is plotted against percentage of 
volatilization in figure 6 
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7. Rates of Volatilization of Hydrofluoro- 
carbon Polymers in a Vacuum 


In view of the complex nature of the thermal 
degradation of polyvinyl fluoride, polyvinylidene 
fluoride, and polytrifluoroethylene, involving split- 
ting off of HF, in addition to scissions of the chain, 
activation energies would be of little significance. 
However, the shape of the rate curves might be of 
interest as revealing details of the mechanism of 
thermal degradation. One rate curve for each of 
these polymers is shown in figure 7. The important 
aspect of these curves is not their relative rate, for 
this is shown in figure 2, but their shape. The poly- 
vinvl-fluoride curve beyond 19 percent of volatiliza- 
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of thermal de gradation of 


a function of entage 


tion resembles that of a first-order reaction, wh 

the polvtrifluoroethylene curve resembles that of 
zero-order reaction, at least in the range 25 to 8 
percent of volatilization. The rate curve for pol 
vinylidene fluoride is conditioned primarily by t! 
stabilization effect, above 50 percent of volatilizatior 


8. Discussion of Results 
8.1. Teflon 


In considering the mechanism of thermal degrad; 
tion of Teflon, the following experimental facts shoul 
be borne in mind: The monomer yield is almost 
percent; the materia: softens and slumps above abo 
50 percent of volatilization; the reaction is of fi 
order; and the rate of volatilization is very likely) 
dependent of chain length, because Teflon, tet: 
fluoroethylene photopolymers, and air-heated Tefk 
all had the same rates of volatilization and the sar 
volatile products 

On the basis of these facts, a mechanism involving 
unzipping of monomer units at free-radical ends « 
chains is assumed. The absence of fraction IT wou 
seem to indicate that the kinetic chain length is ver 
long, so that once unzipping is initiated, it continues 
in most cases to the end of the molecule. Initiatio 
may take place when free radicals form either throug 
breaking off of foreign elements, or groups of el 
ments, at the chain ends, or when a break occurs 
the chain due to thermal agitation. Such thern 
breaks are more likely to occur in the long chau 
than in the short chains, and result in free-radi 
chain ends. The general trend during pyrolysis 
for the average chain length to become short« 
indicated by the fact that the residue softens at abo 
50 percent volatilization. On the other hand, t! 
shorter the chains, the greater will be the tenden 
for their recombination at their free-radical ends 


8.2. Hydrofluorocarbon polymers 


The substitution of one or more hydrogen atom: 
for fluorine on the chain changes radically the n 
of the polymer. Unlike polytetrafluoroethylene, non 
of the hydrofluorocarbon polymers yield any appr 
ciable amount of monomer. Instead of monome! 
the volatiles consist of HF and chain fragments ©! 
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ssizes. The following mechanism is suggested 9. References 
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An Analogue Computer for the Solution of the Radio 
Refractive-Index Equation 
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1. Introductory Concepts 


wue and digital computation devices are not 
ghtening the load of mathematicians and 
ngineers, but also are in many cases enabling 
arch worker to explore problems that have, 
vast, been neglected because of the volume of 
tions required prior to embarking on the 
s of the actual problem. Such a case exists 
-propagation studies in the form of atmos- 
fractive-index computations 


Development of the Basic Refractive-Index 
Equation 


wo parts of the equation for the radio re- 

index as developed by Debve [1]? are: (a) a 
roportional to density resulting from the 
d polarization of the gas molecules in an 
al field; and (b) a term proportional to the 
ent dipole moment of the molecules 

atest microwave determinations of the con- 


of this equation yield the following ex- 
- 


77.6 


1S 10¢ 
T (p+ 7 ); 


refractivity of air, 

refractive index of ai 

absolute temperature in degrees Kelvin 
total pressure in millibars 

partial pressure of water vapor in millibars 


[he Refractive Index in Radio-Wave-Propaga- 
tion Studies 


angles of elevation, the curvature of a 
to a high degree of approximation, 
» the gradient or change of refractive index 
ght. Consequently, the radio horizon fof 
HF communications is determined solely by 
active-index gradient. With an in 
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Laboratory 
and has superseded other method 


the refractive-index gradient, the distance to the 
radio horizon frequently can increase fourfold and 
under extreme conditions may become theoretically 
infinite Therefore, a knowledge of the refractive- 
index profile with height would suffice for refrac- 
tion studies If, however concerned with 
reflection from elevated layers, then one needs the 
additional knowledge of the distribution of refractive 
index throughout the layer 

In radio surveying and navigation 
absolute value of the refractive index is 
of the fact that the velocity of electro- 
magnetic waves in air is inversely proportional to 
the refractive index. 

For a typical propagation study, at least 
Weather Bureau radiosonde ? stations are selected, 
one on each side of the radio-transmission path 
Each radiosonde records data at between 15 and 20 
different elevations. For 1 year’s radio data, there 
would be at least 21,600 calculations of N to be 
made for the two radiosonde stations, taking obser 
vations twice each day. Therefore, it can be seen 
that a great deal of time is consumed in laborious ceal- 
culations prior to making a radio-propagation study 


one 1s 


the 


needed 


studies 


because 


two 


1.3. Present Methods of Calculating the Refractive 
Index 


One of the most common means of reducing the 
calculating time is the utilization of nomograms. <A 
recent example of the nomogram approach is that 
of Swingle [3], which utilizes a separate nomogram 
for each term in the refractive-index equation and 
has a maximum possible accuracy of +1 N unit 
The that utilize a single nomogram are 
never more precise than Swingle’s method, and it is 
not uncommon to encounter errors of +5 N units 

Another approach to solving the equation is to 
use specially constructed slide rules. Weintraub [4] 
utilizes 2 slide rules, 1 for each term in the refractive- 
index equation, and attains an accuracy of » N unit. 


solutions 


The radiosonde is & measurin 
ture, pressure, and relative | 
level to 100,000 ft above sea level 


umidity of 








A major factor to be considered is that both the 
nomogram and the slide-rule methods require the 
striet attention of the calculator, and, as such, fa- 
ticue and the human element restrict not only the 
average accuracy of the results, but also the number 
of consecutive computations possible at one time 
The time necessary .tO compute the separate terms 
n the refractive index and the subsequent addition 


l 
using slide rules, 


varies from approximately 
to from 30 to 45 sec, using nomograms 

The material to follow describes a relatively sim- 
ple electric analogue computer that has been designed 


30 sec 


to solve this equation 


2. Computer Requirements 


The general requirements that a computer designed 
to solve this equation must meet are much the same 
as those that would be placed on any computer. An 
ideal approach to the problem would be to make the 
computer as accurate, as fast, and as convenient as 
However, it would be well to take a more 


possible 
limits to 


realistic approach toward finding practical 
each of these three requirements 


2.1. Accuracy 


In general, the computer accuracy need be no greate! 
than is consistent with the accuracy of the available 
data to be used. Striving for greater accuracy would 
only lead to greater expense and perhaps undue com- 
plications in the circuits. However, the possibility 
that the accuracy of the available data may at some 
future date be increased must be taken into account. 

At present, the greater part of the available 
meteorological data has been obtained from radio- 
sonde ascents. The accuracies associated with these 
data would dictate a computer accuracy of no greater 
than 1 or 2 N units of refractivitv. However, future 
plans indicate that some of the data will be obtained 
from relatively accurate recording instruments placed 
at discrete elevations on specially located steel towers. 
To accurately solve the refractive-index equation, 
using data obtained in this manner, it will be neces- 
sary to specify an accuracy of + N unit 


2.2. Speed 


A very obvious limitation on the overall time 
required to complete one computation is that the 
total time can be no less than the time required to 
supply information to the computer, which is propor- 
tional to the variables in the equation. This is 
commonly referred to as the “setup” time. In 
operation, the operator will be required to read data 
sheets, which have temperature, pressure, and rela- 
tive humidity recorded as numerical values and 
introduce this information into the computer. Using 
either of the two more common media of dials or 
pushbuttons, it would appear well to allow at least 
3 sec for each of these operations. Allowing another 
3 sec to complete the solution would indicate a figure 


of approximately 12 sec to obtain one solutio 
the computer This figure, if met by the con 
would be a considerable improvement over thx 
obtained by either the nomogram or the slid 


2.3. Convenience 


In designing a computer of this type, whicl 
perform a large number of similar computatio: 
much cannot be laid on its convenier 
operation. A machine that has been designe: 
this in mind not only speeds up the operationa 
for the operator, but also greatly reduces opera 
errors. Certain general features emerge from a 
ful consideration of the convenience requirement 

The dials or pushbuttons used to introduc 
variables should be arranged in the same rel; 
order as the data on the data Also 
should be convenient to read and should be calibr; 
in the same units as the original data. The m: 
of reading the answer should be such that the oper 
tor may read and record the result~from the san 
convenient position as used to set up the comput: 

With these general requirements in mind, let 
next investigate the basic electric circuits that w 
be required in an electric analogue computer design 
to solve the refractive-index equation. 


stress 


sheets 


3. Basic Computation Circuits 


If the constants 77.6 and 4810 in eq. (1) are | 
placed by K, and A, respectively, and e by its equiv 
lent, e,#7, the saturated vapor pressure times t! 
relative humidity, the equation will now appear as 


N = (p+K, 7 RH): 


It is to be noted that for the present e,/7' will | 
considered as a variable, as are 7, p, and RH. | 
discussion to follow describes the method used | 
eliminate ¢,/7 as an independent variable and thu 
reduce to three the number of variables that mus 
be introduced into the computer. 

Inspection of the equation reveals that multiplica- 
tion must be accomplished in obtaining the product 
Kz (e,/T)RH. Addition is next used to add p to 
this product, and finally, division is utilized to divid 
this binomial by 7/K,. The circuits to follow ind- 
cate how each of these three mathematical operations 
may be duplicated electrically. In each of thes 
circuits e, represents a d-c voltage source, with zer 
internal impedance, and é is the output voltage, # 
measured with an infinite impedance voltmeter. Th 
variable-resistance symbols represent linear potenti- 
ometers so designed that the voltage as measured 3! 
the slider is equal to the percentage rotation times 
the voltage applied to the potentiometer. X is used 
to represent the percentage of total rotation *% 
measured in the direction of the small arrow ass0- 
ciated with X. 
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3.1. Multiplication 


ircuit shown in figure 1 is presented as a 
| of multiplying two variables by a constant 
uation that the ratio of output 
© input voltage in this circuit is 


expresses 


X,RR 
X,f.+R, 


Ri X,)R.+- 


the equation ® may simplified as 
R sie 
——— XoX 

R,+R,° 
juation shows that when /??,/(?,+-R,) is made 
» the desired constant and X, and X; equal 
variables, the voltage ratio, ¢9/e;, will equal 
duct of the constant and the two variables. 


3.2. Addition 


method of adding a constant, a variable, and 
stant times a variable is shown below. R,, 
are identical potentiometers mounted on 
ie shaft. Their sliders are so phased that 
oth are on zero at the same time. X, may 
ised to represent the percentage of rotation 
potentiometer. The subscripts chosen for 
sistance values will prove to be convenient 
when combining the basic circuits into the 
bined computation circuit. 
tio of output to input voltage is given by 


R, R, 


> , Rea tX Ra 
RAR, Ast ' 


T R,) P R, 


R,-+ R:) (Ra 1 Rs) 
Ry R, T | T Rs 


and R, R, ; | ? Re’ 


the solution obt sined by e« +) st ws 
This maximum error occurs at a value of 


$s 20 times PR, 
f 0.07 percent 


4 Maximum 
2 equal 


or, rewriting, 


Se 


ia) > T 
Rh, LR 


If the correct values are chosen for the resistances, 
this circuit (fig. 2), as shown by eq (4 
to solve the addition of a constant, a variable, and 
a constant times a variable. It is to be noted that 
the measured at the input terminals 
is a constant and equal to PR, 


, Can be used 


resistance as 


3.3. Division 


A method must next be found to divide a constant 
by a constant plus a variable. The circuit shown in 
figure 3 accomplishes this electrically. Again, sub- 
scripts have been chosen as a convenience to later 
work 

The equation for the voltage ratio is 

R, 

?; X R 


R,/R; 
+ R;)/Rs+-X 


Equation (5) is of the proper form to solve a divi- 
sion problem of the type encountered in the solution 
of the refractive-index equation, that is, a constant 
divided by a constant plus a variable 
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Ficure 3. 








4. A Circuit to Solve the Refractive-Index 
Equation 


4.1. Limits on the Variables 


Before proceeding to combine the three basic cir- 
euits. it will be well to investigate each variable 
separately in the equation and discover the limits 
that may be placed on each of them In addition, 
the variables in the equation must be modified so 
that they may be interpreted as the Y variables as 
the This requires that the 
maximum value of each referred to a 
value of unity, and as such can be made equal to a 
percentage of rotation as expressed by the Y vari- 
Consideration will also be given toward estab- 


used in basic circuits 


variable be 


ables 
ishing limits on the variables such that standard 
dials graduated (or marked) with 10 major divisions 
to the revolution may be used 

The temperature, 7, as found in the 
represents the absolute temperature and 
written as 273+-¢, where ¢ is in degrees centigrade 
Because the data obtained from the recording instru- 
ments are in terms of degrees centigrade, it will be 
more convenient to have the dial marked in these 
same units. Investigation of the data currently 
available disclosed that temperatures to be encoun- 
tered will lie between +50° C. This range permits 
each revolution of a 10-turn potentiometer to rep- 
resent a change of 10 deg, and the dial can be marked 
to read zero at its midscale position. With regard 
to the circuit, the temperature expression may be 
rewritten 7'=100(2.23+-X, As a check on the 
validity of this expression, consider the absolute 
temperature at 0° C. X,, or the fractional rotation, 
would equal 0.5, and the expression would give an 


temperature of 273° K, as was to be 


equation, 
may be 


absolute 
expected 
Initial investigation of the pressure variable indi- 
cated that a range of 0 to 1,000 millibars would be 
adequate to accommodate all the current and future 
data. However, closer checks revealed that a range 
of 100 to 1,100 millibars would be more suitable 
This range still permits assigning 100 millibars to 
each revolution of a 10-turn potentiometer, but re- 
quires that a constant be added to the circuit to 
take account of the 100-millibar shift. The expres- 
sion for p is written p= 1,000(0.1+-X, Here again 
Y, is the percentage of dial rotation, with the dial 
calibrated from 100 to 1,100 millibars. As an ex- 
ample, with the dial set for 600 millibars, VY, would 
equal 0.5, and p would equal 600 
The relative-humidity variable lends itself ad- 
mirably to expression as a percentage rotation 
Because it varies from 0 to 100 percent in the same 
manner as the X variable, the dial may be grad- 
uated 0 to 100 and used directly. Investigation of 
the data disclosed that the complete range 0 to 100 
would be required to accommodate all the data 
Although the ¢,/7 variable will later be eliminated 
as an independent variable, its limits will be con- 


sidered at this time. Using the 1951 Smithsonian 
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Meteorological Tables [5], ¢./7’ was found to 
from 0.0003 at 50° C to 0.3820 at +50° ¢ 
order to use all the available rotation of a 

potentiometer, it will be necessary to expres 
variable as 0.3820 Xe,/7T and use a specially m 
dial to position the slider of the potentiomete: 


4.2. The Basic Circuits Combined 
With the above information in mind, eq 


be rewritten 


K : s 
, LOOO(O.1 4 ) + 0.3820. ‘ 
100(2.23-+.X,) 10" \ 


or 


, LOK, P = 
N=9.934X; [0.1 +X, + 0.3820 


hk, x 
1000 


The basic circuits will now be combined, and it will | 
shown that, by properly selecting the resistan 
values, the combined circuit can be used to solve t] 
refractive-index equation in the form of eq (6 

Using the equations already established for +) 
basic circuits, the equation for the combined cireu 
of figure 4 may now be written. The output voltag 
of the division circuit is applied as the input voltag 
to the addition circuit, and the variable X, is ir 
cluded in the last term of the equation for the ad 
tion circuit. Justification for this action is found | 
comparing the multiplication circuit of figure 1 w 
the right-hand portion of the addition cireuit 
figure 2 


The equation for the ratio of voltages in figure 4 


R,, 
R Rew , R,R: >: 
Rh. : y ent A Tr RR, T PR.) X2X ] 








Combined computation circuit. 








will show the sim- 
It is now only necessary 
to the 
The output to input voltage 


omparison of eq (6) and (7 
of the two equations 
sate the resistance combinations in eq (7 
ints in eq (6 

s then equal to NV 


limination of e,/T7 as an Independent 


Variable 
T has 


the 
shows 


( onsidered as i 
The curve 


to this point ¢ been 
ite variable in 

in figure 5,* that e./7 has explicit 
lence on the temperatur This fact can be 
» eliminate e,/ 7 as an independent variable and 
juently reduce the setup time for the computer. 
discussed in section 4.1., if a linear potenti- 
r is to be used to introduce the e,/7 ratio, its 
nust have a special marking The same ratio 
be established by using a linear dial and a 
tiometer with an output function that matches 
rve in figure 5. Unfortunately, the curve 

be duplicated by any convenient mathe- 
However, it is possible to ap- 


equation 


expression 
ate the curve by a series of straight lines con 
the circled Although 


proach does not vield a perfect solution at all 


points on figure 5 
the accuracy can be made as high as required 
rely increasing the number of straight-lin 
mations 

ethod known as 
was used to altel 


loading, or ‘‘pad- 
the normally linear output 
)-turn potentiometer This linear relationship 
n resistance and shaft rotation can be altered 


resistance 


POTENTIOMETER ROTATION , DEGREES 
720 1080 1440 i800 2160 2520 2860 3240 3600 
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-10 '¢) 10 
TEMPERATURE , °C 


er 


by connecting padding resistors between tap points 
on the potentiometer coil so that the output will vary 
in a approximation to the desired function 
By properly choosing the resistance of the pads, the 
straight lines venerated between the tap points can 
points 


close 


be made to cross the desired curve at two 
The number of tap that at 
no point would the straight lines deviate 
from the theoretical enough to contribute over }.N 
unit error to the final result 
Examination of the curve 
the following 
the desired 


points was chosen so 


SeTLeS ol 


that 
positions would 
accuracy: SLO, 1260, 1530 
2340. 2610. 2880. 3150. and 3420 \ 
potentiometel! selected for this 
as 2,000 ohms was the lowest resistance, 


alise losed taps 


located at degree 
provide for 
LSOO, 2070 
2? 000-ohm 
modification 
commercially available, 10-turn potentiometer that 
could be obtained for this Approximat 
resistance values were calculated to pad the 2,000- 
potentiomete! to the desired points on the 
and it was found possible to hold the total 
after padding to 200 ohms. The exact 
values of the padding resistances were determined 
bv the voltage-comparison method described below 
4 500-ohm 10-turn potentiometer was placed in 
series with the tapped potentiometer \ voltage 
was applied to the and the 
output voltage at the slider was read, potentiometer 
voltmeter. The ratio of slider volt- 
age was calculated for each tap point, and, starting 
from the lowest tap point the values of the padding 
resistors were adjusted to make the voltage ratio 
check exactly with the value of e./7T. The 
500-ohm potentiomete! after every 
change in padding resistance so as to keep the total 


was 


purpose 


ohm 
curve 
resistance 


series combination, 


ty pe pret sion 


desired 
was adjusted 
series resistance constant and equal to 523.5 ohms 
Because the maximum value of e,/7 was to be 
0.3820. the 500-ohm potentiometet! replaced 
by a fixed resistance (/?;) of 323.5 ohms upon com 
pletion of the padding operation 

This method provided much gvreatel 
than could have been obtained by calculating the 
values for the padding resistors With the type of 
precision voltmeter used, it was possible to read the 


was 


aHcecuracy 


voltages to five-figure accuracy, and consequently 
changed to permit adjusting 
the voltage ratios at the tap points to an accuracy 
ot 0.01 percent Wire-wound resistors of the 
temperature coefficient that used in the 
potentiometer were specially constructed for this 


purpose 


the resistors could be 
same 


wire as 


Using the potentiometer just described, it is now 


possible to graduate its dial from 50° C to 50° C 
in exactly the same manner as the temperature dial 


eliminated and the 
directly to 


In fact, the dial can now be 
e,/T potentiometer can be 
the shaft of the temperature This 
has in effect eliminated the « 1 that it 
now depends directly on the temperature dial for its 


setting. 


connec ted 
potentiomete I 
variable, in 





6. Final Design Considerations Both sides of this equation must be multiplied hy 

1,000 to permit solving for values of N up to | 009 

Investigation of the actual data disclosed th al 

the values of N will be between 0 and 1,000. Th 

As shown by eq (7), the computation circuit & | dia] attached to R, must therefore be graduated fron 

developed to this point will provide a solution to the | © to 1.000. Comparison of eq (6) and (9) requires 
refractive-index equation in the form of a voltage that the following identities be satisfied in or 
A more convenient and accurate method of utilize the final circuit to solve the refractive 


6.1. Choice of Solution Presentation 


ratio 


presenting the solution 1s to combine the composite equation: 

circuit of figure 4 into a bridge circuit and use a P 

» | nce ia( Py 
ala ce. 1000 bagg\ fi 


sensitive null indicator to detect the bridge Ryo 10k, =776 
Figure 6 shows the complete computer circuit dia- ri= 440" 
gram. Potentiometer 2, and the fixed resistance 
PR, constitute one arm of the bridge, and the com- R,+R;) 
hined computation circuit of figure 4 forms the other R 
arm. i, is the same type 10-turn linear poten- 
tiometer as used to introduce the variables. A R 
linearly calibrated dial attached to 2, was used to = 
read N directly. Ry is a wire-wound resistor, as 
are all the other fixed resistances K 

The null indicator used is of a special design - “ : \2 
having a maximum sensitivity of \-in. motion per riz) 1000 
microampere at about the zero position. Its response 
falls off logarithmically on either side of zero and (2, +R) (R. +R 
thus provides automatic protection against excessive ~~ Rit+h+RetR, 
currents caused by large unbalances in the bridge ‘ 


eireuit R. 4 Tt R,. Tt Rea. 


RR, 


There exist several combinations of the al 


6.2. Selection of Resistance Values - : 
resistances that will satisfy the above identities 
The variable resistances or potentiometers must firs 
be chosen on the basis of commercially availab 
written values. In general, they were chosen as small a 
VR possible in order to keep the over-all internal resist- 
. 4 ance of the circuit low with regard to the interna 
Ry+ hi impedance of the null indicator. Low circuit resist- 
ance is necessary in this case to maintain sufficier 
expression may be equated to eq 7 sensitivity of the null indicator. —_ 
R,, and Ry were chosen as 100 ohms each. This 
R(Ro+-R value, by eq (12), required that R,,=10 ohms, whi 


RR, le , RR y X: | (q) | Was still of a size that could be conveniently woun 
%)*~ 
\ 


wy) 


Referring to figure 6, the following expression for 
the voltage as measured at the slider of Ry may be 


With the circuit as shown at null balance, this 


R+R. R.* Xi RAR R, and R, have already been determined as 323 

R ‘ and 200 ohms, respec tive ly. These values aad tl 
value for R,, when substituted in eq (13) give a val 
for R, of 481.00 ohms. Substituting this value 
eq (14), along with the values for R,, R., and R 
gives a value of 5,825 ohms for Ri. Solving eq (15 
gives a value for R, of 591.00 ohms. Equation 
was next investigated in order to determine a val A com 
for R;. As R; cannot have a negative value, J sults o 
must be at least 278 ohms. The closest commer- HH shows th 
cially available value of 500 ohms was chosen for 2 This sat 
Using this value for R:. eq (11) was solved, and f 
was found to be 524.0 ohms. 

Equation (10) was next used to determine the rels- 

tive magnitude of R, and Ry. Solving this equatio 
established the relationship 


R, . 


RR 0.25773. 
‘ 10 


+ 


R,+Ry was chosen to be approximately 2,000 ohms 


Ficure 6. Computer- so that excessive current would not be drawn fro 
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oltage source and yet not be so vreat as to cause 
is loss of sensitivity at the null indicator. This 
for R,+Ry would indicate a value for R, of 
515 ohms A resistance of 500 ohms was 
n for R, because resistors of that value are com- 
ally available. This choice for R, required that 
sual 1.440 ohms. #&, was selected as 50,000 ohms 
ler to satisfy the previously established relation- 
hat R, be at least 200 times R 
ire 6 shows a resistance #,, in series with the 
y supply voltage. This resistor was used to 
e the bridge voltage to approximately 50 volts 
rher voltages a slight heating effect, along with 
sponding resistance changes, was noted in some 


resistors 


6.3. Wet- and Dry-Term Switches 


refractive-index equation may be considered 
ne composed of two parts, a wet term equal 
] KLRHe, T) and a dry term equal to (K,/T p 
of the applications that require the solution 
refractive-index equation also require the 
for either the wet or dry contribution to N 
witches were incorporated into the final circuit, 
permit the computer to solve for either con- 
on and thereby greatly increase the usefulness 
computer 
cure 6, Ri has been divided to form Rg and R, 
s been made equal to 10 ohms and R, equal to 
ohms, thereby keeping the total resistance the 
With the switch connection between Ry and 
sed, the pressure term is effectively reduced to 
ind the computer, after rebalancing, reads the 
rm. Ina similar manner, closing the dry-term 
effectively makes the wet term zero, and the 
ter solves for the dry term only 


7. Results and Discussion 


compute! herein described was built under the 
s supervision at the National Bureau of Stand- 
Boulder, Colo., and is at present being actively 
DY their Radio Meteorology Group The use 
computer has resulted in considerable savings 
over previously used methods 


7.1. Accuracy 


omplete test was run on the computer, and the 
of this test are given in table 1. The table 
s that the error is never greater than + 4 N unit 
satisfies the original requirement If greater 
y were to be required, a greater number of 
oints on the e,/7 potentiometer would permit 
approximation to the desired curve and con- 
tly increase the final accuracy. The need for 
accuracy would also necessitate the 
vith finer graduations, 


use of 


7.2. Speed 


The average operational time of the computer and 
operator was found to be approximately 10 see, 
which is less than the original requirement of 12 se¢ 
This time naturally with the individual 
operator. The shortest time of 8 sec was recorded 
for an experienced operator and the longest time of 


Varies 


13 see 


‘ was recorded for a relatively inexperienced 
operator 


This time could be decreased somewhat 
by using an automatic balance svstem to rotate the 
V dial and thereby balance the bridge circuit more 
rapidly than can be done manually. This refinement 
was not deemed practical for the present computer 
due to the relatively longer time required to set in 
the variables 


7.3. Convenience 


All of the original requirements on convenience 
have been incorporated in the final design. Figure 7 
shows the computer. It is housed in a small sloping 


Ficuri 





front cabinet that can be placed on a desk or table 8. References 
and used in the same manner as a standard calculat- 
Jocause of heating effects, 1t was found P. Debye, Polar molecules, sec. 3.4 (Chemical ( 
Co., New York, N. Y., 1929 
K. Smith, Jr., and S. Weintraub, The constant 
, t : itrs onstan 
bridge circuit continuously With this in mind, a equation for atmospheric refractive index at ra 
foot switch was incorporated in the final design so quencies, J. Research NBS 50, 39 (January 
R P2385 
D. M. Swingle, Nomograms for the computation of 
spheric refractive index, Proc. IRE [3] 41, 385 
and record data 1953 
In figure 7 the various dials are set up to solve for | S. Weintraub, Slide rule for the computation of tI 
620 refractive index of air, Electronics [1] 26, 182 (J 
1953 
Smithsonian meteorological tables, 6th revised 
2 (Smithsonian Inst., Washington, D. ¢ 


ing machine 
desirable not to leave the batter connected to the 9) } 


that the operator could energize the circuit as needed 
and still have both hands free to operale the dials 


V,. with the temperature, 32.0° C; pressure, 
millibars; and RH, 62.5 percent V is read as 277.0 


9c — 
oo oo 


BouLper. Co.o.. Julv 13. 1953 








On the Accuracy of the Numerical Solution of 
the Dirichlet Problem by Finite Differences 
J. L. Walsh’ and David Young * 


This paper derives numerical bounds for the error, in certain closed regions, of the differ 
ence analog of the Dirichlet problem It is concerned only with the difference between the 
exact solution of the difference equation and the solution of the Dirichlet problem. The 
error bounds obtained involve quantities which can actually be computed, such as the mesh 
size, and the oscillation and modulus of continuity of the given function on the boundary 
So far as the method is concerned, the chief novelty is the use of the differences analogs of 
harmonic measure and the Schwarz Alternating Process 


1. Introduction 


Although finite difference methods afford a powerful tool for obtaining numerical solutions 
of partial differential equations, little is known about the accuracy. It is the purpose of this 
paper to derive numerical bounds for the error, in certain closed regions, of the difference 
analog of the Dirichlet problem. We shall be concerned only with the difference between the 
exact solution of the difference equation and the solution of the Dirichlet problem. The error 
bounds that we obtain involve quantities that can actually be computed such as the mesh 
size, and the oscillation and modulus of continuity of the given function on the boundary. So 
far as method is concerned, the chief novelty is the use of the difference analogs of harmonic 
measure and the Schwarz Alternating Process 

Gerschgorin [4] ¢ derived error bounds for boundary value problems associated with elliptic 
partial differential equations. These, and similar bounds derived by Collatz [2] and Mikeladze 
10] involve bounds, in the closed region, of certain partial derivatives of the solution of the 
differential equation. However, the solution of the differential equation itself is not known, 
to say nothing of its derivatives, (although approximate values for the derivatives may some- 
times be found by examining the corresponding difference quotients). Also, it may happen 
that although the derivatives in question are not bounded in the closed region, the solution of 
the difference equation may still converge to the solution of the differential equation 

Rosenbloom [13] presented an error bound for the Dirichlet problem which is closely 
related to Gerschgorin’s but which utilizes special properties of harmonic functions. By use 
of well-known inequalities giving bounds for the partial derivatives of a harmonic function at 
an interior point in terms of the oscillation on the boundary and the distance from the bound- 
iry, he obtains upper bounds for the derivatives in closed subregions. Then, by solving the 
difference equation on subregions and, as the mesh size approaches zero, letting these sub- 
regions approach the given region, Rosenbloom obtains an error bound involving w*(4), the 
modulus of continuity in the closed region of the solution of the Dirichlet problem. However, 
Rosenbloom does not discuss the question of finding w*(4) in terms of w(4), the known modulus 


of continuity of the given function on the boundary. Furthermore, in practical numerical 
work one would not wish to change the boundary of the network as the mesh size is decreased 

In section 3 we use the explicit solution of the difference analog of the Dirichlet problem 
for the rectangle, obtained by Le Roux [9] and by Phillips and Wiener [12], to derive an error 
bound involving bounds for the derivatives of the given function on the boundary. Wasow 
15] has obtained error bounds for the rectangle which are applicable when the boundary values 


Pr Presented to the American Mathematical Society, September 1951; abstract published in Bul. Am. Math. Soc $7, 478 (1951 The paper 
was prepared in part while Professor Walsh was engaged as consultant on a National Bureau of Standards contract with the University of Cali- 
rnia at Los Angeles. The research was also supported by the Office of Naval Research under Contracts N5ori-07634 and N5ori-76, Project 22, 
with Harvard University, and later by the Office of Ordnance Research, U. 8. Army, under Contract DA-36 034-ORD 966 

? Harvard University; consultant, National Bureau of Standards contract with the University of California at Los Angele 

3 University of Maryland 

* Figures in brackets indicate the literature references at the end of this paper. 
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have bounded third derivatives. Our results are somewhat more general in that we allow the 
second and third derivatives to fail to exist at a finite number of points. Also, we obtain error 
hounds for the case where only the second derivative is bounded, and we allow the first and 
second derivatives to fail to exist at a finite number of points. In sections 4 and 5 harmonic 


measure and its finite difference analog, discrete harmonic measure, are used to obtain bounds 


for the moduli of continuity of harmonic and discrete harmonic functions, respectively, in 


certain regions in terms of (6). Phillips and Wiener [12] showed the existence of such bounds, 
whereas we obtain suitable bounds in a precise numerical form. Then, in section 6 these upper 
bounds are used to vield a uniform error estimate for the closed rectangle in terms of the oscillia- 


tion on the boundary and 6(w In section 7 the extension to regions made up of two or more 


overlapping rectangles is discussed 


2. Discrete Harmonic Functions 


Let A and & be arbitrary positive numbers, and let L{A,k] denote the set of points (z, y) 
such that both z/h and y/k are integers. Let 2 be a simply connected closed region with interior 
R such that the boundary S of @ consists of straight lines, each of which is parallel to a co- 


ordinate axis and contains a point of Lih,k]. Let Q, denote the subset of points of L{h,k] 
contained in 2. Two points (2,,y,) and (22,y.) of L{h,k] are adjacent if 


[ (xy; —22) /hAP?+[(yi— ya) /AP=1. (2.1) 


A point of Q, is an interior point of Q, if the four adjacent points belong to Q,. All other points 
of 2, are boundary points. We let R, and S, denote, respectively, the set of interior and 
boundary points of 2,. Evidently, we have R;,SR and S,CS. 


A function U(z,y) defined on Q, is said to be diserete harmonic,’ (d. h.) in R, if it satisfies 
the difference equation 
&[U (x,y)]=[207/(1+ 0) |[U(r+h,y) + U(e—h,y) —2U(z,y)] 

+-(2/(1+ 07) |[U(a,y +h) + U(a,y—k) —2U (z,y)]=0 
where 
o=k/h. 

The finite difference analog of the Dirichlet problem is the following problem: given a 
function f(2,y) defined on S, to find a function U(z,y) defined on Q,, d.h. in R; and coinciding 
with f(z,y) on S;. The existence and uniqueness of a solution of this problem for bounded 
regions is easy to prove; see for instance Gerschgorin [4]. The convergence to the solution 
of the Dirichlet problem has been proved using nonconstructive methods by Le Roux [9], 
Phillips and Wiener [12], and others 


3. Error Estimate for the Rectangle Under Differentiability Assumptions 


Let 2 be bounded by the lines z=0, r=a, y=0, and y=6, where a= Ah, b= Bk and A and 
B are positive integers. Let f(x,y) be defined and continuous on S and let u(z,y) and U’(z,y) 
denote respectively the solution of the Dirichlet problem and its finite difference analog with 
boundary values determined by f(z,y). In this section we shall derive an upper bound for the 
error U(z,y)—u(z,y) in the region Q,, under certain assumptions about the derivatives of f(x,y). 
If f(x,y) is defined and continuous on S, we denote generically by u;(z,y) and U;(z,y) the solu- 
tions respectively of the Dirichlet problem and its finite difference analog with boundary values 
f(x,y). 

First, one can verify directly that the function u,(z,y) defined by 
u,(z,y) =f(0,0) +[f(a,0) —f(0,0)] (z/a) +-[f(0,b) —f(0,0)] (y/b) 

+[f(a,b) +f(0,0) —f(a,0) —f(0,b)](z/a)(y/b), (3.1) 


* Heilbronn [5] introduced the term discrae harmonic function, and studied the properties of these functions 
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which is linear on every line parallel to either coordinate axis, is both harmonic in R and d.h. in 
R,. Thus if f,(z,y) is defined on S and equal to u;(z,y) on S, then u, (x,y) = U,(2,y) in Q, 
Next, if fr(z,y)=f(z,y)—f(z,y) on S, then we have by linearity of Laplace’s equation and 


Evidently both U,(z,y) and u,(z,y) vanish at all the corners of @ 

We study arbitrary boundary values f(x,y) by studying in turn four functions each of 
which has the boundary values /,(z,y) on one side of 2 and values identically zero on the other 
three sides. By the symmetry of the situation, it is sufficient to study in detail only one of the 
latter functions. Let us set /;(z,y)=f,(z,0) when y=0 and f,(x,y)=0 elsewhere on S. Since 
f,(z,y) is continuous and vanishes at the corners of Q, f;(z,y) is continuous on S 

It can be verified directly ° that if the Fourier series of f, r.0) converges to fa(2,0), then we 
have . 

sinh \(n3b a)(l y h 


sin (nrz/a ; 
sinh (nrb/a 


sin (nrt/a)dt, 


and that 
sinh (mab a)(l y 4 


sin (nrz/a 
sinh mab a 


i 
At=(2h/a) >} sin (nwjh/a)f,(jh,0), 


— 


and where m and n satisfy the relation 


sinh (mrk/2a asin (nrh/2a 


Now let g(a f.(x.0 We now prove the following theorem 

THEorEM 3.1 Let g(x) satisfy the following conditions: (a) g(x) and its first (s—2) deriva- 
tives are continuous, O<a2<a;: (b 7 D(x) fails to exist or fails to be continuous for, at most, a 
finite number of points in the interval 0<2<a, and g“~"(x) is uniformly continuous in each open 
interval in which it exists and is continuous: (c g(r exists, erce pt possibly for a finite number of 
points, and is bounded, O<2x<a 

Then un iformly for all L,Y) Alp, we have 


OC)D(h/a 
if s=2, and 
Ce) Dy(h/a)’, 
3, where 


1/24) x?(1+ 0") coth (arb/a)/ae, 


D,=2(a/xr)*[2J,M,/a+M,], D,=2(a/r)*[2(J.+ 1)M,/a4+- Msg] 


l 


2/xc) sinh (o 


a \<é 
—— 


* See, for instance [8], pages 95 and 96 
’ See, for instance, Le Roux [9] or Phillips and Wiener [12 





and J. 4 BRS denotes the n umber of points in the open interval O<2r<.a at which g (x) does not 
1.2.3. denotes the least upper bound of the modulus of g (x) 


s not continuous Tlere V/ 


erst or 7 


PROOF 


LeMMa 3.1 lf m and n satisfy (3.7), then 


2 ee 


o~ renth a 


PROOF Let p ) q=nrh/2a. Then (3.7) takes the form 
sinh p=ce sin q. (3.14 


We now study the function p(q) defined by (3.14) in the interval O<qSr 2 One can 
verify directly that p(0)=0, p’(0)=e, and 
a(1-+o*) sin q 
o* sin’ q . 
— o°)(1—2e* sin*g) cos g 
i) ( 7 - : = , 
f |+o* sin*qg)""* 
where primes denote differentiation with respect to q By the extended mean value theorem 


Therefore, m<n. On the other hand, we also have 


since p”’ q) * o. we have pseq 
_— a*)| 1—2e* sin*q 
p vie ; 
o* sin*g 
since |(1—2o*sin’g) (1 + o¢°sin“q 
Avain using the extended mean value theorem, we get 


1/6 gp’ g ((Q). 


and 


Pq oq S (1/6 qo 1+ o 
and the lemma follows 
LEMMA 3.2 lf m and n satisty 


Ps and it a is give v by (3 12), then 


a 


Proor. Defining p and q by (3.14), we have 


d . 
7, (P(9)/a)=|P'(g)a— P(q)l/ 9°. 

q 
By the mean value theorem, we have p(q)=qp’(&) (0<&<q). Since, as shown in the proof 
of lemma 3.1, p’’(qg)S0, we have p’(q)<p’(&) and p’(q)q—p(q) <0. Therefore, the ratio 
p(q)/q is a nonincreasing function of gq, (0<q<-/2). Its minimum value in the interval 
when qg=x/2. Therefore, p(q)/q2 p(x/2)/(#/2)=(2/x)sinh~'(¢), and 


O<qSr/2 is assumed 


the lemma follows 
Now for convenience, let us define the function 
sinh [(mxb/a)(1 — y/d)| » is 
, (3.16) 


sinh (mrb/a) 


We shall study its behavior as a function of both m and y, where m is assumed to be a con- 


tinuous variable 
LEMMA 3.3 fi m= 0, then for 0< y* bh we have 


d 


00> 
: dm 


l.(Y¥) > —(ry/a) coth (mrb/a) exp [—mry/a]. 
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Proor. Differentiating (3.16) with respect to m we obtain 


y b cosh [(mrb/a)y l —y¥ b)] sinh mab a 


rb a 
1—y/b sinh~*(mxb/a 


— cosh mab a sinh [ mab a 


rb al ly y/b) coth [(mxrb/a\i— y/b6)|—coth (mxrb/a 


We note that for K>0, z>0 

S (9K + 9k 

coth (Kr Kr esch? (Ka 1 sinh 2K x ; 2K 3 
2 sinh*( Kz 


l ; 
(x coth Ka 
ly 


But since sinh(2Az)>2Kz, the last expression is nonnegative; therefore by the mean value 


theorem we have 


y b coth| mrb/a)(1 y b |<coth mab/a 





On the other hand, 


cosh|(mx6/a) (1 y b ro ve 
coth(mrb/« 


d } 
y TO yi0 
: sinh(mrb/a 


Dn | 
dm _— 
cosh[(mab/a) (1 y/b)) 
rt a coth mrb a) y b ae Loe 
( cosh(mxb/a 


sinh (mry/a 
rh a coth m rh a y b ; Y Vi b Pn (Y) 
: sinh (m2b6/a) cosh (mrb/a) . tel, 


the brackets are, respectively, nonnegative and nonpositive, whence 


For 0<y<6, the terms in 


ry/a) coth (mxrb/a) T,,(y). 


Finally, we note that 
(l—y b } 


exp [—2 mab/a) 
Sexp (—-Mry/da), 


- y exp mary a 
2 1—exp |— 2mrb/a} 


provided® 0 <y<b The lemma now follows. 
3.4. lf m and n satisfy (3.7). then 


LEMMA 3 
d ; 
| ls ”| <S coth arb a ae | n 
dm , m | 


verified that if 0<?t, O<A, then {te 
The lemma follows from lemmas 3.2 and 3.3, and from the fact 


where a is given by (3.12 
Proor. It is easily Ktl < (Ke Therefore, 
y exp(—may/a) <[(mr/a)e]~. 
that coth(z) is a decreasing function of z, for > 0. 
Lemma 3.5. If m and n satisfy (3.7), then 
Pn (y l.(y)| <C(h/a)*n’?, 


where C and a are given by (3.10) and (3.12), respectively. 


* This inequality was proved by Phillips and Wiener [12) 
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Proor. By the mean value theorem we have 


> d 
. — |" . ‘(2 - m—T 
I ™ Y I n\l Ma x m<m,<n E : d 1. my 


The lemma follows from (3.13) and (3.17 


Lemma 3.6. If A, is defined by (3.4), then 


A,| <D,/n', 
where D, is defined by (3.11 

Proor. The proof involves the use of repeated integration by parts of the integral ex- 
pression for A,, 


and is similar to that given by Jackson[8, p. 13-14] for the case s=2. We 
omit the details 
We now define the function /,(z,y 


i 
féz,y)=0 unless y=0 and f,(7,0)=30 A, sin (nrz/a 
n=\ 

For s=2,3 and 0<z<a, we have 


en 
A, sin 
i 


naz/a)i\<D, >\ n“<D, t—*dt—D 
n=A+l 


f,(x,0) —g(x)| <[D,/(s—1)] (A/a). 
On the other hand, if we replace f,(z,0) by f,(z,0) im (3.4) and (3.6), 
1 ,2,...,4-1 Therefore, for (x,y) « L {h,k], we have, since sin (Arz/a) 
U (z,y)—ulz,y ’ A, sin (nwxz/a)(l,,(y)—T, (y)], 


U (z,y)—u,(z,y)| < >> \A, 


— 


A-1 
‘Tn(Y)—T.(y)| <CD,(h/a)’ >) 

n=l 
by lemma 3.5 and (3.18). It can be shown that 


n* 


A h 1 
atc) ' 


~ -log A 1 —log A/a 
Therefore, if s 


Uy (x,y u,(r,y)| <CD,(h/a) 


(3.23 
Evidently (3.22) could be used to obtain an error bound for the case s=3, and this was done 
in the original manuscript of the present paper. 


Since then, however, a paper by Wasow [15] 
has appeared which contains a stronger result for the case s=3 than the original form of theorem 
3.1. A combination of his methods and the original ones now yields an improvement (the 


present theorem 3.1) of Wasow’s result in this case since we allow g(z) and g®(z) to fail to 
exist for a finite number of points 


Thus following [15], we obtain from lemmas 3.1, 3.2, 3.3, and the mean value theorem 


Ply l.(y)| <Cae(h/a)?(xy/a)n®exp(—anry/a). 
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By (3.20) and lemma 3.6 we have 


U,(x,y t4(2,y) <CDyae (h/a)* (xy a) >) exp anry/a)<CeD,(h/a 


since for y>0,>> exp annry/a jexp (ary/a 1]-'<(ary/a 


Since d.h. functions possess mean value properties, they assume their maximum and mini- 
mum values on the boundary as do harmonic functions. Since the functions [us(2,y) —u4(2,y)] 
and [U’;(2,y) — U’,(z,y)] are harmonic and d.h., respectively, we have 


Max |q (a f.(x,0)!, s(x, "(z,y)| <Max |g (a 


U,(z,y 


+2 Max g(a 


3.24), the theorem follows 
4. Harmonic Measure 


Let 2 denote a simply connected region with interior 2 and boundary S. Let S’ denote a 
subset of S consisting of a finite number of connected subsets of S. We define the harmonic 
measure, (h.m.), H[(z,y),S’,Q] as the unique function which is harmonic and bounded in R&, is 
continuous in 2 except perhaps at a finite number of points of S, and equals unity on S’ and 
zero on S—S’. The properties of harmonic measure have been studied in considerable detail, 
see Nevanlinna [11], chapter II] 

By analogy we define discrete harmonic measure, (d.h.m.), for regions of the type described 
in section 2, as follows: H,[(z,y),S’,Q;} is a function d.h. and bounded in R,, equal to unity on 
SNS’ and to zero on (S—S’)QS, 

For bounded regions the existence and uniqueness of h.m. is well known, see [11]. The 


existence and uniqueness for the half plane can be proved by the use of conformal mapping 

The existence and uniqueness of d.h.m. for bounded regions follows from the existence 
and uniqueness of the solution of the difference analog of the Dirichlet problem. Later we shall 
prove existence and uniqueness of d.h.m. for a half plane and for certain other unbounded regions. 

We list some elementary properties of d.h.m., which are analogs of well-known properties 
of h.m. 

0<H,[(z,y),S’,Q,] <1 for all S’SS. 
If VX SS”’CS, then 
H,|(z,y),S’ 22) < AI (z,y),S8”’’,Q,). 
" 


If S’ is included in the boundary of both Q,; and QF where Q,SOF and if (z,y)eQ,, then 


H,| L,Y S’,2,) < Hi [(z,y) ,S’ QF . 


The first property follows at once from the fact that the maximum and minimum values 
of d.h. functions are assumed on the boundary. The second follows from this fact and from 
the fact that the expression { 7,[(z,y),S’’,Q,]—Hz|(z,y),S’,Q:]; is nonnegative on S,;. The 
third follows since by the maximum and minimum principles H;{(z,y),S’,Q7]20 and hence 
H,|(2,y),S’ ,Q7|—H_[(z,y),S’,Q:)}2 0 for any point (z,y) on the boundary of Q,. This is the 
so-called principle of gebietserweiterung. 


We shall use h.m. and d.h.m. for two purposes. In this section and in section 5 lower 


bounds will be derived to enable us to obtain upper bounds for the modulus of continuity of 
harmonic and d.h. functions in a closed region in terms of their moduli of continuity on the 
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In section 7, upper bounds will be derived which will enable us to use the Schwarz 


boundar 
Alternating Process and its difference analog for overlapping rectangles. 
Let 2 denote the rectangle of section 3, and let S’ denote the side contained in the line y=0. 


THEOREM 4.1 If (2-a/2 y?< 8, then 


v(b/a) (6/a 
v(b/a) (6/a 
v(b/a) =[4+ 2’coth?(xb/a)]'”. (4.6 


PROOF Let u(z,y), U(2,y) be harmonic and d.h. functions respectively in ? vanishing 


on S—S’ and equal to sin(rz/a) on S’. By (3.3) to (3.7) we have 


Ss b sinh[{(mxb/a)(1— y/6)) 


sinh[(4b/a / ’ 
and »! Ssin( rsa 


a 
sinh(rb/a sinh(mrb/a 


where sinh(mk/2a) =e sin(rh/2a 
By lemma 3.1 we have m‘ Also by lemma 3.3 we have (d/dm)T,,(y) <0. Therefore, 
y)>T,(y), and we get 
>u(r,y). (4.7 


Now, for 0<@<* » have® sin @2>(2/r)@; and hence: for 0O<2<a/2, we have 


| r But since sin [(r/a)(a/2+p)]|=sin [(x/a)(a/2—p)] for all p 


sin(wr/a)7 2r/a 


it follows that sin(rr/a)> 2—r\, (O0<rSa 


Also 


sinh{(9rb/a)\(1—y/b r/a)cosh{(#b/a)(1 — y/b)}| <(/a)cosh(x b/a 


dy 
and by the mean value theorem 


sinh|[ (xb/a y/b)| >sinh(rb/a ry/a)cosh(xb/a 
Therefore, 


l—u(zr,y) : ry/a)coth(xb/a)] } < (2/a)\a/2—z\ + (ry/a)coth(rb/a 


By the Schw arz inequality the last expression does not exceed 
y*] 2y(b/a < bla)vlb/a . 
provided I a/jZ) y 


We now observe that //[(z,y),S’,Q|>u(z,y), Ay[(z,y),S’,Q,)> U(z,y). The theorem now 
follows from (4.7 

Now let 2; denote the subset of LjA,k] such that y>0 and let R, and S, denote the interior 
and boundary of Q,, respectively 

Tueorem 4.2. For the region Q,, d.h.m. exists and is unique 


PRooF Consider the sequence Hy LY)}, where 
HH (z,y) =H, {(z,y) S’OSY? Qe], 


and where 2°’ denotes the rectangle with vertices (nh,0), (—nh,0), (nhnk) and (—nh,nk). 
Evidently 7’ (z,y) exists, since d.h.m. exists for bounded regions. By (4.3) it follows that the 
sequence { //;"’(z,y)} is nondecreasing and bounded above by unity. Therefore, a unique limit 
exists which we denote by H,(z,y 


* This is Jordan's Inequality; see Copson [3, p. 136 





Since 77" (x,y) is d.h., it satisfies (2.2 Taking limits of both sides of (2.2), we find that 
H,(2,y) satisfies (2.2) and is therefore d.-h. Evidently H/, rey) <1, Ay(a.y 1 on S’OS, and 
H,(zr,y)=0 on S’—S, Hence //,(2,y) has all the properties required for d.hom 

In order to establish uniqueness, we prove 

LEMMA 4.1 Tf Ula y) is d.h., bounded in R, and vanishes on S,, then U(xr.y) vanishes in R, 

PROOF By hypothesis there « xists a constant P such that for (z y é82,, we have |U'(a I P 
Given any (#,y)eQ2; and any €>0, let a and 6 be chosen so that (a,b)eZ{h.k) and 


Let 2* denote the rectangle with vertices at (ab ' a,b), (a.0). and 1.0). and let 7, denot« 


the interval a<sr<a Since lU(z,y) <P on S* and vanishes on /,, we hav 


By theorem 4.1 we have, replacing a by 


HH, 


provided r } . ) But by hypothesis, we hav , U Substituting we 


H 


1/2): hence 


[his proves the lemma. 

The uniqueness can now be proved by assuming two d.h.m.’s and showing that theit 
difference vanishes identically 

The proof of theorem 4.2 is complete 

We note that the existence ol a unique bounded solution of the difference analog of the 
Dirichlet problem for 2, with bounded boundary values is almost immediat: If the boundar\ 


values are determimed by g(x), then the limit of the absolutely and uniformly convergent series 
>» H, 4 ph.Q lg ul 3s d.h. in R,. bounded in 42, and equals giz) ons Thus a solution 


xists. The uniqueness follows at once from lemma 4.1 
We next consider the region z20, y20. To find d.h.m. for subsets of the line y=0 we 
perform a sign-changing reflection about the line z=0 and use theorem 4.2. Since the d.hom 


is zero on the line z=0, this can be done in such a way that the new function will be d-h. in 


R Similarly to find d.h.m. for subsets of the line y=0 we reflect in the line z=0 


By similar methods the existence and uniqueness of d h m. can be established for othe 


regions suc h as the semi-infinite strip Usrsa, y > () 


5. Modulus of Continuity 
The modulus rf continuity of a function / roy) in & closed region {2 1s define dl by 
LUB\ f(ay,y:) —f (xo¥: 


where (4 Yi) (22,Ye) & and (2 I2)° y Y2)* Sh Evidently if 6 is not less than the diamete 
of 2 then w(4) equals the oscillation of f(2,y) in Q 

For harmonic functions we prove the following theorem 

THFOREM 5.1 Let 2 denote a bounded simply connected closed reguor with interior R and 
whose boundary S is a closed Jordan curve with the following property: there exist constants 
ro>0 and @>0 such that for any point P of S there exists a circular sector with vertex at P, with 
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radius ro and included angle @ containing no point of R. Letu r,y) be harmonic in R. continuous 


nQ with modulus of continuity w(6) on S If for some positive ha to eee 
Max!ry.D.D (ro/D) tt!" 
then the modulus of continuity of u(a.y) in Q satisfies the inequality 
w | D(6/D)*!* +] + (4. M/x) (6/D)*! "+9 


where ¥ Max 27 


Max |u rey Nin [u(r y 


’ yes 


N the Ot lat 7] 
PROOF 
LEMMA 5 ; deno ‘Ses, ¥ and let C, denote the region 


> 
7] p J ()} 


PRoot It is easy to show that if ¢ is the angle at r,0) from J, to the cirele C through 
the points r.0), (r,0) and (z,y), then H7[(z,y),J,,.y2 0] =1—¢/x 
Now if p/r 1—cos ¢)/sin ¢, then (, is contained in the region bounded by C and the line 


0. Hence, by the minimum principle for harmonic functions we have for (x,y) e( 


FH (x,y) ,,, , 2/r)tan , 2p/ rr, 


and the lemma is proved 


y- <r, O< tan yt) <W <2 and let B, , denote 


LemMa 5.2. Let C,y dei gion 2+) 
the hound 1 rad ) r ' cr then B J denotes the line 0<2x< , y=0 Ii par and if 


i 


Proor. We first consider the case Y=. If 2°+y7*’=r* and y>0, then as in lemma 5.1 


we have Hi (z.y)./..y>0]=1/2 Hence, as we verify at once, 
H{ (x,y), 1,,C,)=2H 1 (z,y),1,,y > 0] 
If (nye, -. then by lemma 5.1 we have 


H 


The lemma can now be verified for the general case by mapping (,, onto C,, by means of the 
conformal transformation (w/) */* We omit the details 

Now by the maximum principle for harmonic functions, for given 6 the maximum value of 

y U'(2,ye)! for (2,—2 Yi — Yo)? <*& occurs when either (x,y) or (2%,y2) belongs to S 
If (2,,y)eS, then for all r such that 0<r<ry there exists a circular sector C,, containing at 
least one point of R with vertex at (2,,y,) with radius r>6 and angle y. such that B, , the union 
of the bounding radii is disjoint from FR. Since the theorem is trivially true if (22,y2)eS, we 
assume (2%,y2)e2. Let Q, denote the closure of the connected component of (22.4%) for the region 
QNC,y. Evidently QeC,,. Since for (z,y)e2,NS, we have 


H r,y),B, y,C,. 9] Hi (x,y 20S8,2,] 


.¥ 


and for (2,y) contained in Q,; and on the are of C, 4, we have 





then for all (2,y7)eQ, it follows that 
Fl (x, 


But by lemma 5.2 we have 


HH| f 2, 


since the lemma is obviously true if we rotate the sector 
Therefore, we have 


w(r) 1 (22,42) ,C,,40S,2)] + MHA] (22,42) 8 7,4, 02] 4.\f/x) (6/r)*/9 


Now by the assumptions on 6 we have D(é/D)**® <D(ry/D Hence, we can choose 
D(6/D)*/"*® <ro, and the theorem follows 
We now consider the case where Q is an arbitrary simply connected region. We first prove 
the following general theorem, which is essentially a formulation for harmonic functions which 
s equivalent to Carathéodory’s form of the theorem of Milloux for analytic functions 
THEOREM 5.2 Let G be a Jordan subregion of |z\<1, whose boundary consists of a Jordan 
ire a, which passes through 0 plus an are a, of T:\z\=1. Let u(z) be harmonic and bounded in 
G. continuous in the corre sponding closed region except at the end points of a, and a, equal to unity 
the interior points of a, and to zero in the interior points of a» Then in ¢ ery po nt 2of G we hare 


> | t/r)tan 


Proor. Let be conjugate to u(z) in G, and let / exp|—u |. Except 
perhaps at the end points of the ares, on a, we have | f(z e~', and on a we have |f l; on 
these open ares |/ is continuous in the two-dimensional sense. Then by the form of Mil- 
oux’s theorem presented by Carathéodory, [1], 354, we have for zeG 

log! ‘ + (4/7) tan 
is we were to prove 

THEOREM 5.3 Let 2 denote a bounded simply connected closed region with interior R and 
houndary S Let u L,Y be harmonic in R. continuous in 2 and having modulus of continu ty 

6) on S lf Dis any positive constant, then the modulus of continuity of u(x y) in Q satishes the 
nequality 
w[D(6/D)*?] 4M 3) (6/D 


w all 6<D. 


Proor. As in the proof of theorem 5.1, it is sufficient to show that if C, is a circle 


with center (2,,y,) on S and radius r and if (22,y2.)eQ2 and (2;—.2, Yi — Yo)? <& <r’, then 
FT| (22,Ys 


The theorem then follows if we let r=D(6/D)'* since 6<D and hence D(é/D > 6 

Now, if 2 is a Jordan region, (5.4) follows at once from theorem 5.2, since tan 
lO! <1 We indicate the modifications necessary to include the case where Q is an arbitrary 
simply connected region 

Let O be a boundary point of Q, let C, be the circle r, and let 2) be a point interior to Q 
and to ¢ If no point of © lies exterior to C,, then a function v(z) harmonic in ? and continuous 
in 2 equal to unity on S is identically unity in 2 so we have v(z i. We-proceed to study the 
contrary case. If points of P lie exterior to C,, such points can be joined to 2 by a Jordan are 
lying wholly in RF, so at least one are of C, lies in R. Let the totality of mutually disjoint arcs 
of C, in R be A), Az, Denote by /, the subregion of ? interior to C, containing 2) and 
by 2, the sum of Jo, its reflection (inverse) in C,, and the ares A, which form part of the boundary 
of Ry; if an are of C, is part of the boundary of PF, it does not belong to PR We modify 2, by 
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adjoining to FR, the interior of each Jordan curve that can be drawn in #2; and to which O and 

Denote this new region by 2. Then R, contains 2, is simply connected, and 
point. The subregion 2, of R, interior to C, is also 
of C, of which every point is 


are exterior 
has O as either an exterior or a boundary 
simply connected, and part of its boundary is an open are A 


accessible: we choose Ao as the largest such are, so that every boundary point of 2; not on Ag is 
or is a point of C, not contained in an are of C 


either a boundary point of FR interior to C, 


consisting wholly of accessible boundary points of PR 
in which Ay necessarily corresponds to an arc of y) it is clear that 7(2,Ao,R 


By a conformal map of 2; onto th 


interior of a circle > 
exists and is unique this function takes the boundary value unity at every point of A» and th 
boundary value zero at every boundary point of R not on the closure of A Carathéodory s 


proof of Milloux’s theorem is valid for the region FR, the distance from 2 to the boundary of 


RP. is not greater than so (5.4) is valid for an arbitrary simply connected region, and th 


theorem follows 


Che modulus of continuity of a function F(z,y) on any subset Q, of L[h,é] is defined by 
Li B f r,.V 


where ry ro. é\. and where (4 J y Yo)" 0 
For d.h. functions in a rectangle we prove the following theorem 


Let {) denotl the rectangle O<7<,.a, US y< A. where ad eL\h,k] Let 4? 


THEOREM 5.4 
lenote the subset of L h j conta ed nm If { ry is d h an ie. the interior of $2, and has 


winuity «@ 0 A ’ boundary of Qz, then for 6<7 ab the modulus of 


nN odulus 


nequality 
hk] M/2)v(1/2 


where V a iy } ( } f ».2 and 4.6), respectively 
Proor. <As in the proof of theorem 5.1 the maximum value of |U’(x,,y ¥2,Y2) for 


I } / y.)? <& is assumed when either (x,,y;) or (2,y%2) belongs to S, Let us assume 
that (2,.y,)eS, Now if (22,¥2) also belongs to S; the theorem is trivial, since 6<, We there- 


fore assume (22,y.)eR 


and in luding one of the sides of the rectangle, 


Let p denote a straight line containing (7,,y 
») 


Let C denote the closed interior of a semicircle with center at (x ¥,), with radius 6 Z ov 


h+k}, with bounding diameter included in p and containing at least one point of R,. Since 
5<r we have 6,226. Evidently there exists a rectangle 7 included in C with one side con- 
tained in p and with vertices contained in L{h,k] such that the sides perpendicular to p are 
at least half as long as thos« parallel tO p and such that the latter sides have length at least 
>( O7 


By (4.2) and (4.3) we have 


Now, for (2,y)eS 
r y 
Moreover, since the intersection of p and the boundary of 7/NQ, is contained in TNS, we have 


H,{(a,y), TO p,T|=0 for all points of the boundary of 2,9 7 not contained in S,NT. Therefore, 
for all points of the boundary of 7NQ, we have 
Hil (a2, S,,.Q,0 7) = H,\(a,y), Tp, 7). 
It follows that for all (x,y)e7NQ, we have 
H,\(2,y), T0S2,2:0 T)=> A, (2,y), TN p,7). 
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Now since (22.% we have by theorem 4.1 


i 
i 
U (x,y 


d the theorem follows 


In the above proot i / A’, one could have obtained a res tangle with the des ead properties 


tting ©’ have radius 2 h) Therefore, we hav: 


(‘OROLLARY Tt h ke the 


W/2): 


6. Error Estimate for the Rectangle In Terms of the Modulus of Continuity on 
the Boundary 


In this section we obtain an erro! bound for the re tangle of section 3 under the assump- 
that the function f(2.y which determines the boundary values has modulus of con- 


ons The function f,(a y defined in se tion 3 is also continuous, and we denote 


tv w/fd 
its modulus of continuity for 0 << 


Let / ry) and ua(r.y) be given by 3) and .5 


a, y=0 
respectively We define the function 


by the partial sum 
> a, 
the . ire the Fourier coefficients for f.(c.0 g(x) as defined by $4) and where the 


ents d 


lote respec tively the solutions of the 
If / 0, the values are determined by f.(a We now prove 


are summation coefficients defined by Jackson [{7, p. 9]. Let w;(2,y) and U(x. 


Dirichlet problem and its difference analog vanishing 


S except fo! y) () 


THEOREM 6.1 If (2.yeQ, and 0 <2: <y<hb. ther 


SA*OX h a ea 
VW*=Max 


where ¢ and a are 


PROOF It can be 


nrxz/ayr 


narra, y 


and n» are related by }.7), and where ls Y) is defined by 16 Evidently we have 


6.6 


By lemma 3.3 we have 








By lemma 3.2 we have 7 ri; <A, and 


ry/a) coth (arb/a) exp anmy/a). 
Using lemma prool of lemma 3.5 
Cae(h/a)*(xy/a)n® exp anry/a 


Therefore, we 


' « 
Ss‘ j : ; <Cae(h/a)(ry/a > n exp (— anry/a), 


which is convergent for y>0 


Now the following statement can be verified easily: if G(z) is a continuous nonnegative 


funetion for c>0. nondecreasing for 0<2<.Y, and nonincreasing for X <z, then 


5G n)< Gitjdt+GXx 


Cherefore. since for vy >0 we have 
Max [z° 
it follows that 


S*n® exp ann f r yl ( 3a/amrye)’ <6(a/ary)* 3a/ am Ye 


Therefore, we a 


AC(h/ay y/a 


For the Fourier coefficients, since |g(t)|<.M*, we have 
i(t) sin (nat/a) dt 


Jackson 6] proved that 


From 6.6), (6.9), the theorem how follows 


COROLLARY . , Q, ( Sra, eSySb we have 
8. *( » hia)? é/d 2Kw 
where A wa posui / onstant less than 


PRooF. 
{ ry t3(\7,Y))* ai) UY 
By the maximum and minimum principles for harmonic and d.h. functions, we have for (z,y)«Q, 
U’s(z,y) — Us(2,y)|,\us(z,y) —us(z,y)| < Max! f,(x) —g(r 
0<s<e 


Using a theorem of Jackson [7, p. 7], we obtain 


Max | f, (a g(x) < Kw 2a/A), 





here A isa positive absol unt less than Che corollary now follows irom 
Wi shall how cderive 


he boundary valu 


theorem 
bound for the original probli m whet 
COoTTe! 


or al mig the 


THEOREM 6.2 





a appre 


PRroo! 


LEMMA 6.1 Tf u,(a2.y lefined by (3.1), then for all (x,y) we have \u(x.y 
Proo! since wry and 4,(2,y are 


ssumed on S We can assume 


int whose ordinate is zero 


harmonic, the maximum of 


loss of generality that the 
By (3.1) we have 


u(r.y 
maximum occurs 


without 


riqajnulog 


u(a,0)| << Max u(z,0 Min u(z,0) <M, 


3 


nd the lemma follows. 


The above result cannot be improved as one can show by considering the case where 
ry) is identically zero and u(z,y)20 
From lemma 6.1 we conclude that |us(z,y)|< M. 


In particular we have |g(z)|< M/; hence, 


Mt<M. 


6.20 
LeMMA 6.2 


6M da. 6.21) 
Proor. Evidently we have a 


where w,(6) denotes the modulus of con- 





tinuitv of u,(2.0) considered as a function of z. But since 


r.0 6/a)iu(a,o u“(OO)], 


we have w,(6)<6\M/a, and the lemma follows. 


Lemma 6.3. If e<r<a—e,e<y: @. 


16M XC, (h/a ( (A : + w(2k)+-2.M(h/a+k/b)] 


7 rr . r . 
Proot We observe that r.) x u(x,Yy The function [U.(2,y 


sum of [U’5(z2,y i,(x,y)| and three other terms of similar type. The lemma 


/ ry Is the 
follows from (6.12), from (6.20) and from lemma 6.2 We note that from the hypothesis 
wi have ed.eé h< 


Lemma 6.3 affords us an error bound for those points of 2, which are at a distance of not 


less than ¢« from the boundary In order to obtain a uniform error bound we now consider those 
points which are within ¢ of the boundary 

Since 4/r<v(1/2)/2 it follows from theorem 5.1, (with D=r), and theorem 5.4 that the 
moduli of continuity of U'(2,y) and u(z,y) do not exceed the upper bound for w7(6) given by 
theorem 5.4, provided ab If (x,y)eQ, and if (x,y) is within ¢ of a point of S, then ther 
exists a point of S; at a distance not greater than ¢e from (z,y We remark that every point of 


Q, is within (ab)'? of some point of S Since U'(x,y)=u(z,y) on S; we have 
Qwe (€ 
If / denotes the right membe f 2), then we have 
Max T ;2wz € 
In order that the error in §2, should approach zero with the highest power of + 
whenever w(5)/5 is bounded as a function of 6, we choose e=r(hk/r?)?”’. Evidently we have 
‘(h/a)*"(b/a V7 (hk /b)*7* (ald 
The theorem now follows from (6.23), lemma 6.3, and theorem 5.4 


For the special case of the unit square theorem 6.2 and theorem 5.4, corollary, give 


COROLLARY If a=b=—a=1, then we hare 


u(z,y)|< Max! J;;J2}, 


whe re 
») 


32.M\Ch?" + 8K[w(2h 2Mh\, '=(x*/12ae) coth (ar) ~.572, 
2w[2'/2(h?/? h)] My(1/2)h?", 2 m)sinh '(1)~.561, 


27 ») 
,~3.566 


Substituting numerical values, we obtain 
J, <66.Mh2" + 240(2h) +48Mh, J. <2w[2"/* (h?’7 + h)]+4Mh2", 


The above expressions for J; and J, represent a slight improvement in the formulas previously 
given by the authors [14] 


7. Other Regions 


In this section we consider regions of the type described in section 2 other than rectangles. 
The case of two overlapping rectangles is studied in detail as an illustration of the method 
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which can be extended to more complicated regions. We first prove a theorem valid for two 
overlapping regions whether rectangles or not 
THEOREM 7.1 Let Q2=0'U0""’, where R'NR’ is not empty, and where neither region 
neludes the other Let u(a y) and Ua u denote, re specti ely, the solution of the Dirichlet proble m 
nd its finite difference analog for Q Let 2,, 2. and 2 denote the points of Lih Jd helongina to 


’, and Q’’, re spectively Let ir Sy R; ri oy R; / 
Max H, 
riyeT 
Nezt. let lL’ and U"’ de note functions d 0; rie 0’ yo ai to ( SY and S 
espectively Then for (2 ye 


< Max (A,/ 


Max |/ 


y)eR 


Max / 


r.y)eR; 
PROOF Let 


Max |U LY) UT, Y 


T.yleT; 
Clearly, for LLY "yp 
y 


uw, Max |U (2,1 
r.y)eT 
Hence, by the maximum and minimum principles for d.h. functions F”’ 
Similarly, F’’ Sus” E’ Therefore, 


By the maximum and minimum principles, 


Max |U(z,y 
Ty) ef; 
and the theorem follows 
The proof of this theorem was, of course, motivated by consideration of the difference 
analog of the Schwarz Alternating Process. Indeed, in order to solve the difference equation 
for Q;, one might guess values for U’(z,y) on S,’N/,’’ and then solve the difference equation 
in R,’ obtaining in particular values on S,’’NR,’. Next one solves the difference equation 
for R,’’, using the computed values on S,’’NRP,’, and obtains new values for S;’NR,’’. This 
process can be repeated and the successive values thus obtained converge to the exact solution 
of the difference equation. Moreover, the rapidity of convergence can be estimated if 4, and 
uw. are known. In fact, one can show that after a complete iteration the maximum error is 
reduced by a factor of (1—n,y. 
If for 2’ and 2” the quantities Z’ and E”’ are known whenever the modulus of continuity 
and oscillation of u(z,y) are known on S’ and S’’, then by using theorem 5.1 or theorem 5.3, 
E’ and E”’ can be computed, provided the modulus of continuity of u(z,y) on S is known. If, 


moreover, u, and/or uw. are known, an error bound for the composite region can be obtained. 
We consider now the case of two overlapping rectangles 2, and Q, with sides parallel to the 
coordinate axes. If the intersection of the interiors is not empty, and if neither region is a 
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other. then the cases A through F shown in figure 1 occur 


lhe composite regions of cases C and D are of the same type as those of cases A and B 


respectively. The problem of determining 4; and yw, become problems in d.h.m. for the rectangk 


ol the types shown In figure 


























FIGURE 2 


Upper bounds for d.h.m. are required in points on the closed dotted lines for the open ares 
on the boundary indicated by heavy lines 


By (4.3), these upper bounds are not greater than the upper bounds for the problems show 
in figure 3 
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PROBLEM I PROBLEM IT PROBLEM IT PROBLEM IV 
(FOR TYPE A) (FOR TYPE 86) (FOR TYPES B, ANO F) (FOR TYPE E) 


FIGURE 3 


For problem I, since the harmonic measure of the open arc is required, }is an upper bound 
on the dotted line by symmetry 

THeoreM 7.2. Let Q; be bounded in part by a segment of a line l. Let R, contain no point 
of l Let l’ denote a line perpend cular to l and contain ing a point of INS, Let S’ de note any 


subset of INS, contained in one of the open half planes bounded by l’. If (2,y) is any point of 
RAV, ther 


H, z,y),S’,Q,) < 


Proor. The theorem follows easily by (4.2), (4.3), and problem I. 
As a corollary we have an upper bound of 3 for problem II. Also, for problem III if we 
extend the lines whose harmonic measure is required to the right, we get by symmetry as it 


360 





problem I } as an upper bound for the d.hom. on the left dotted lin Similarly 


is ah luppe r 
bound for the d.h.om. on the right dotted line 


For problem IV the situation is somewhat mors complicated Wi 
THeoreM 7.3: Let Q denote the region: 2>0, y>0, and let Q, 
ntained in Q Let I,. Ly denote respect ely the intervals 0 

t (aA) be longs to Lihf We have 


prove 


/ I 


FIGURE 4 








Here Yo the smallest number not less than a/2a such that y kis an intege) The quar ities a and ea 
are defined by (3.12) and (2.3 respectively. 


a 


PROOF Let J*, J* denote respectively the intervals x >a, y=—0, and y >h a 0) We first 
prove 


LL.eEMMA 7.1 r O<y. then 


PROOF By theore m7 
HN (x,b).y 


Hence, we certainly have 


y bNz 0 } y < . >()) 


But by a sign-changing reflection about the line y=0, we can show that the left member of the 
above expression minus the right member equals /7/,|(z,b) 77 ,Q,])—Hz[(2,6) 5,21) 2 0 

This proves (7.3 Evidently (7.4) can be proved by the same method 

From lemma 7.1, we conclude that 


HeA(2,b) 1, ,2,) < HI (7,6),7 4,22) 1/2)H,| (2,6), 0< yNx=0,Q,). 








1/20H 
Cherefore, the d.h. funetior 


VW 


satisfies the conditions 


W(x,0 ral 7; 0<y); Wia,y)S1/2, (OSy 


Let ¢Q)* denote the r rial l ( > 0, and let OF denote the corresponding subse 


Lih,k). Evidently 


W here 


We now prove 


Lemma 72 sb ’ nd bounded for ryeRt and such that 


[ 2(O,y U. a,y 0, (y2U 


< Max [w;l w)c/Yo). 


Proof ’ 5 3.6) and taking the limit as 6 becomes infinite we get 
) sin (nrz/a) exp mry/a), 


where Ah a 


By lemma 3.2, we have m= a Also, since cot [nx/2A] <2A/nx for n< A we have 


: l 
‘ , . rT exp Tray/da) 
> > 4/xn) exp ( ‘ > 1/rn) exp (—rayn/a 2/x) log 2 . 


1—exp ray/a 


Evidently if y>a/2a, we have 


Now since 0<w<. 1, the function 
l1—w)Yy/Yo 7.12 


is dA in RZ and is not less than U,(z,y) for SiN y < yo and for y= yo. 


Therefore, for y <y we have 


r,y < ] (1 W)Y/Yo, 7.13 





and the lemma follows. 


Now U, roy 2T rey 


and the theorem now follows 

In applying the theorem to problem LV, we note that we can assume the value zero at the 
pomt (0,0 without changing the d.h.m 

Thus, an error bound can be obtained for the case of two overlapping rectangles. The 
above methods can also be used to find error bounds for other regions of the type considered in 


) 
section 2 
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